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Since quantum field theories do not possess proper position observables, Unruh-DeWitt detector 
models serve as a key theoretical tool for extracting localized spatiotemporal information from 
quantum fields. Most studies have been limited, however, to Unruh-DeWitt (UDW) detectors that 
are coupled linearly to a scalar bosonic field. Here, we investigate UDW detector models that probe 
fermionic as well as bosonic fields through both linear and quadratic couplings. In particular, we 
present a renormalization method that cures persistent divergencies of prior models. We then show 
how perturbative calculations with UDW detectors can be streamlined through the use of extended 
Feynman rules that include localized detector-field interactions. Our findings pave the way for the 
extension of previous studies of the Unruh and Hawking effects with UDW detectors, and provide 
new tools for studies in relativistic quantum information, for example, regarding relativistic quantum 
communication and studies of the entanglement structure of the fermionic vacuum. 


I. INTRODUCTION 

While the spatial and temporal meaning of wave func¬ 
tions in first quantization is straightforward, it is notori¬ 
ously difficult to extract spatiotemporal information from 
the quantum fields of second quantization. This is in large 
part due to the fact that there are no position observables 
in quantum field theory (QFT) [1]. As a consequence, 
one often works with nonlocal field modes instead, for 
example, in order to calculate S'-matrix elements in the 
case of particle physics or Bogolubov transforms in the 
case of QFT on curved space. 

Actual experiments, of course, probe quantum fields on 
finite patches of spacetime. A technique which does allow 
one to model the extraction of localized spatiotemporal 
information from quantum fields was pioneered by Unruh 
and DeWitt [2, 3]. The key idea is to probe quantum fields 
by coupling the quantum field in question to a localized 
first quantized system, called a detector. An excitation of 
the detector system is then interpreted as the absorption 
and therefore detection of a particle from the quantum 
field. 

For example, a hydrogen atom, first quantized, can 
serve as such a detector system for the photons of the 
second quantized electromagnetic field. The technique 
of probing quantum fields by coupling them (at least in 
gedanken experiments) to so-called Unruh-DeWitt (UDW) 
detectors has been used successfully, in particular, to an¬ 
alyze the Hawking and Unruh effects. One of the key 
insights gained from these studies has been the finding 
that and why the very notion of particle is observer de¬ 
pendent. 

More recently, UDW detectors have been used exten¬ 
sively in studies on quantum communication via field 
quanta [4, 5] and, more generally, in studies on a host 
of effects related to the presence of spatially distributed 
entanglement in the quantum field theoretical vacuum 


[6-10]. In spite of these successes, however, most of these 
studies have been limited in scope to the model of a UDW 
detector that probes a quantized real scalar boson field 
through a linear coupling. 

Here, our aim is to widen the applicability of the general 
UDW detector approach. To this end, we investigate 
Unruh-DeWitt detector models that probe bosonic as 
well as fermionic, nonscalar fields through both linear 
and quadratic couplings. We show how certain persistent 
divergencies can be overcome through renormalization 
and we show how perturbative calculations with UDW 
detectors can be streamlined through the use of extended 
Feynman rules that include the scattering of detector 
excitations. 

Historically, idealized particle detector models were 
first introduced by Unruh [2], in an attempt to resolve 
the well-known ambiguity in defining field states corre¬ 
sponding to physical particles on curved spacetimes or for 
noninertial observers [2, 11, 12]. This was the finding that 
different quantization procedures yield incompatible Fock 
spaces and thus lead to ambiguous definitions of parti¬ 
cle number eigenstates [13, 14]. Unruh’s groundbreaking 
idea is best summarized in his now famous dictum “a 
particle is what a particle detector detects.” He initially 
modeled a particle detector for a quantized scalar field 
as being a quantized scalar field itself, albeit one that is 
restricted to a cavity so that its excitation events come 
with spatiotemporal information. Unruh calculated that, 
when uniformly accelerated through the Minkowski vac¬ 
uum, this detector will measure a flux of particles that 
are thermally distributed, which established the famous 
Unruh effect [2, 14, 15]. 

A few years later, DeWitt improved and simplified Un¬ 
ruh’s original model by introducing approximations that 
effectively replace the detecting field with a nonrelativis- 
tic two-level quantum system ]3], inventing what is now 
called the Unruh-DeWitt detector. The detector couples 



2 


directly to the field <I> through its monopole moment /r: 

i7int oc . (1) 

This model fixed minor problems with Unruh’s initial 
suggestion [13]. More importantly, the UDW detector 
has the advantage of being much easier to work with in 
calculations, because the detector is described by quan¬ 
tum mechanics rather than quantum field theory. Other 
detector models for quantized scalar fields followed, us¬ 
ing both two-level systems and full fields as detectors, 
and featuring different ways to couple the detector to 
the fundamental quantum field (see e.g. [13, 16-18]). In 
the present paper, we will generally call any model that 
uses DeWitt’s monopole as the detecting system a (UDW- 
type) detector, irrespective of the field it probes and the 
coupling it uses. 

UDW particle detector models have proven to be very 
versatile and useful tools. They were first used to analyze 
the effects of accelerations, spacetime curvature and hori¬ 
zons on scalar quantum fields (see e.g. [2, 19]). Further, 
UDW detectors have been used in quantum field theory 
to quantify vacuum fluctuations and the structure of vac¬ 
uum states ]20, 21]. In relativistic quantum information, 
pairs of detectors are crucial for measuring and harvesting 
the entanglement of quantum fields ]22-25], and serve as 
senders and recipients of information in quantum signal¬ 
ing [5, 26]. Moreover, they can be good toy models for 
light-matter interaction in quantum optics ]27, 28]. 

Because UDW detector models for quantized scalar 
fields have proven so valuable, the question arises to what 
extent analogous techniques can be developed and applied 
to other types of fields, such as quantized spinor or vector 
fields, which in Nature are of course more common than 
quantized scalar fields. This question is of fundamental 
interest by itself and there are also concrete applications 
where detector models for quantized spinor fields could 
immediately be useful, such as for analyzing the fermionic 
Unruh effect ]29], and for investigating the related Hawk¬ 
ing radiation of Dirac particles [30]. Moreover, it has 
recently been pointed out that such more general UDW 
detector models could resolve ambiguities that arise in 
defining measures for the entanglement between disjoint 
patches of a fermionic field ]31]. 

A key aim of the present work is, therefore, to find a 
counterpart of the UDW detector for quantized spinor 
fields. More precisely, we are interested in a pair of particle 
detector models, one each for quantized scalar and spinor 
fields, that is comparable in the sense that differences in 
the detectors’ reaction to scalar and to spinor fields are due 
to the probed field alone, not caused by any other feature 
of the respective detector models. We will focus on only 
half-integer spin fermions and integer-spin bosons, i.e., we 
will here not be concerned with nontrivial spin/statistic 
combinations that may arise in low dimensional systems. 

To the best of our knowledge, only two detector models 
for quantized spinor fields have been suggested so far: 
the first one by Iyer and Kumar [32] uses a quantized 
scalar field in a cavity as detecting system, close in spirit to 


Unruh’s original suggestion. The second one by Takagi ]20, 
33] seeks to imitate DeWitt’s simpler model by coupling 
an UDW-type two-level system to a quantized spinor field 

through 

Hint oc , (2) 

making it a likely candidate for the closest spinor field 
equivalent to the scalar field UDW detector. While the 
above interaction may not be immediate from first princi¬ 
ples, that is, from the Standard Model of particle physics, 
we may for example think of it as an attempt at a first- 
quantized, simplified version of a second-quantized cavity 
detector: in quantum electrodynamics, a spinor field (elec¬ 
trons) is coupled to a vector field (photons) through 

Hqeb oc . (3) 

In this logic, certain modes of the electromagnetic field 
restricted to a cavity with appropriate boundary condi¬ 
tions (realized, e.g., as superconducting mirrors) could 
serve as a detector for the electron field. The spatiotem- 
poral profile of the detector is then determined by the 
extent of the cavity and the time over which the elec¬ 
tromagnetic field is observed. It is conceivable that by 
a series of approximations—e.g., neglecting all but one 
mode in the cavity and restricting its occupation to at 
most one quantum—we could arrive at a model of the type 
of Eq. (2). A two-level system coupled through Eq. (2) 
may then be thought of as an approximation to a cavity 
detector, much in the same way that the UDW detector 
is a simplification of Unruh’s original cavity setup. In the 
following, however, we will keep with the original idea of 
such detectors as simple tools to obtain spatiotemporally 
resolved information about quantum fields, and set aside 
considerations of how exactly they could arise from first 
principles. This is justified by the fact that, as discussed 
above, UDW detectors provide a very useful tool for ex¬ 
tracting spatiotemporal information from quantum fields, 
at least in gedanken experiments. 

We begin in Sec. II with an overview of the UDW- 
type detector models that we will study. In Sec. Ill, the 
excitation probability for the different detector models at 
rest in Minkowski vacuum is calculated for these models, 
also to test to what extent these models are comparable. 
We demonstrate that any UDW-type detector featuring 
an interaction Hamiltonian which is quadratic in the field 
has divergent excitation probabilities. This is true for 
spinor field models like Eq. (2), but also for quantized (real 
and complex) scalar field detector models with coupling 

Hint oc . (4) 

Unlike divergent probabilities reported before ]13, 20, 34, 
35], it is not possible to regularize these divergencies 
by considering a detector of finite size and switching it 
on and off adiabatically. However, the divergencies are 
readily understood from a field theoretic point of view: In 
Sec. IV we show that they are closely related to tadpoles 
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in quantum electrodynamics, and this analogy leads us 
to a renormalization scheme for quadratically coupled 
detectors. Sec. V is dedicated to the comparison of the 
vacuum excitation probabilities of the detector models. 
We establish that it is indeed justihed to compare the two 
models Eq. (2) (for quantized spinor fields) and Eq. (4) 
(for quantized complex fields), at least at leading order 
in the coupling constant. As an example, their response 
to massless quantum helds in (1,1) dimensions is studied 
explicitly, using both a sudden switching of the detector, 
as well as a Gaussian switching function, in addition to a 
Gaussian spatial detector profile. 

The second part of this work is dedicated to the deriva¬ 
tion of fully fledged Feynman rules for all detector models 
discussed, valid on Minkowski spacetime. The Feynman 
rules are to facilitate investigating whether the discussed 
models are finite, or at least renormalizable, and whether 
they are comparable at higher orders as well. Sec. VI 
summarizes the required calculation methods, such as 
the Feynman propagators and Wick’s theorem, with the 
proofs moved to Appendix G for better readability. After 
applying these methods in Sec. VII to the probability for 
a detector to remain in its ground state when interacting 
with a held in the vacuum state (the vacuum no-response 
probability), the Feynman rules then follow in Sec. VIII. 

In Appendix A, we summarize canonical results in clas¬ 
sical and quantum held theory to introduce notation as 
well as for the convenience of the reader. Appendix B ex¬ 
plains the origin of the different normalization conditions 
for scalar and spinor held mode functions, which plays a 
crucial role in the UV behavior of the theories. 

Throughout this work, natural units c = \ = h are used. 
The metric of (n -I- l)-dimensional Minkowski spacetime 
has signature (-I-, —, • • • , —). 


II. UNRUH-DEWITT-TYPE PARTICLE 
DETECTOR MODELS 

Unruh-DeWitt particle detectors are two-dimensional 
quantum systems. We choose the energy levels to be zero 
and some SI > 0 so that the detector’s Hamiltonian, which 
generates translations in its proper time r, reads 

H |e) (e| = H . (5) 

Here, the ladder operators 

= |e) {9\ , cr" = \g) (e| (6) 

are expressed in term of the ground and excited states 
|g) and |e). We parametrize the detector’s trajectory 
x{t) = (t(r),a:(T)) by its proper time. Since the detec¬ 
tor’s trajectory will be prespecihed rather than dynamical, 
spacetime translational symmetry is broken and overall 
energy and momentum are not conserved (unless one 
considers also the agent that keeps the detector on its 
trajectory). The detector’s effective spatial profile, which 


may be thought of as the distribution of its charge, cen¬ 
tered around x, will be described by a function 

p{x, •) : H” ^ H (7) 

that is normalized to one unit of charge. 



where n is the number of spatial dimensions. For example, 
for a pointlike detector 

p{x,y)=6{x-y) , (9) 

while for a detector with Gaussian profile 

p(a;, j/) = (27rg^)~"/^ ^ . (10) 

It will sometimes be convenient to combine the switching 
function x a-nd the spatial prohle p into a spacetime profile 
/ by defining 

f{x) = f{t,x) = x{t)p{0,x) (11) 

in the detector’s rest frame. 

An UDW detector is to detect held quanta by becoming 
excited. To this end, the held must couple to an operator, 
p, of the UDW detector that does not commute with the 
detector‘s Hamiltonian Without restricting generality, 
one can dehne this so-called monopole moment operator 
of the detector to be 

p = a~. ( 12 ) 

Let us now systematically consider how the monopole 
moment p of an UDW detector can be coupled to the 
various types of helds. 

A. Linear coupling 

We begin by considering a pointlike UDW-type detector 
coupled to a quantized real held 4> through the interaction 
Hamiltonian 

Hi^,=Xp^x), (13) 

where A G R is the coupling strength. This is the particle 
detector model originally proposed by DeWitt [3]. It can 
be viewed as an idealization of the two lowest states of 
an atom interacting with the electromagnetic held. 

As is well known, UV divergencies in the excitation 
probabilities arise in this setup, in particular, if the de¬ 
tector is abruptly switched (i.e. coupled) on or off; see, 
e.g., [13, 34, 35]. The pointlike structure of this detector 
also leads to UV divergencies; see, e.g., [13, 20[. The di¬ 
vergencies can be regularized by switching the detector 
on and off gradually, i.e., by using a sufficiently smooth 
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switching function x(^)) by endowing the detector 
with a finite spatial profile p: 

Hint( t) = Xxir) fj, [ p{x{T),y)^{y)dy . (14) 

JR" 

Let us now consider the coupling of an UDW detector to a 
quantized complex field, <I>. In this case, one cannot couple 
them as in Eq. (13) because <I>^ 7 ^ $ would imply that the 
interaction Hamiltonian is not self-adjoint. In principle, 
one may instead view <& as composed of two quantized 
real fields ® ^ + i‘^’ 2 ) and couple any real linear 

combination of $1 and $2 to the detector. This approach, 
however, would single out a direction in the complex 
plane and would therefore make any C/(l) symmetry for $ 
impossible [20]. The same incompatibility arises with any 
unitary symmetry if $ carries higher-dimensional complex 
group and/or spin representations. We conclude that the 
coupling of UDW detectors to quantized complex fields 
should not be linear. 


B. Quadratic coupling 

Any interaction Hamiltonian between an UDW detector 
and a field must be a self-adjoint scalar. Since an UDW de¬ 
tector’s monopole moment is a self-adjoint Lorentz scalar, 
UDW detectors always need to couple to an expression 
in the field which is also both scalar and self-adjoint. In 
the case of a quantized spinor field, the simplest self- 
adjoint Lorentz scalar is TT = 'I'lq^'I'. This suggests the 
following interaction Hamiltonian: 




real 


complex 


spinor 


- coupling - 

linear quadratic 


1. Hint ~ 4> 

2. Hint ~ 4>" 

- 

3. Hint ~ -1>+-I> 

- 

4. Hint ~ 


Table I. Overview of the particle detector models discussed 
in this paper - listed is the power of the field to which the 
detector (i.e. y) is coupled. 


tool for investigating what difference it makes to an UDW 
detector’s behavior whether it is coupled to a bosonic 
(scalar) or to a fermionic (spinor) held. This is because to 
this end we can now compare model Eq. (16) and model 
Eq. (15) which are both quadratically coupled, i.e., they 
differ essentially only in the spinor structure but not in 
the type of coupling. We will separately investigate the 
impact of quadratic versus linear coupling, such as model 
Eq. (16) versus model Eq. (14). 

In Table I, the four different detector models we will 
discuss are summarized and numbered for easy reference: 
model 1 for quantized real helds couples linearly through 
Eq. (14); model 2 couples quadratically to quantized real 
helds; model 3 to quantized complex; and model 4 to 
quantized spinor helds with the Hamiltonians Eq. (16) 
and Eq. (15), respectively. 


Hintir) = Xxir) y [ p(a7(T), y) T(y)T(y) dy . (15) 

JR" 

Note that the interaction is now quadratic in the held. 
This is plausible also because an UDW detector coupled 
linearly to a fermion held would be able to violate fermion 
number conservation by creating and annihilating individ¬ 
ual fermions. Coupling an UDW detector quadratically to 
the held ensures that the detector can only pair create or 
annihilate fermion antifermion pairs. An example of this 
type of coupling was hrst considered by Takagi [20, 33[ 
and has been used since; see [36-39]. 

It is straightforward to couple UDW detectors quadrat¬ 
ically not only to quantized spinor helds but also to quan¬ 
tized (real and complex) scalar helds while preserving 
symmetries: 

Hintir) = Xx{t) y [ y(a:(T), y) $'f(y)4>(y) dy . (16) 

JR" 

This type of coupling was hrst suggested by Hinton for 
quantized real helds [17]. It may be possible to justify 
interpreting this model as an effective description of some 
fundamental interaction, for example when the scalar held 
describes composite particles at sufficiently low energy. 
However, in the present work, we will use it as a theoretical 


III. FIRST ORDER: EXCITATION 
PROBABILITY IN THE VACUUM 

For simplicity, only UDW-type detectors that are at rest 
in n+1 dimensional Minkowski spacetime will be consid¬ 
ered in our investigation of the different detector models. 
In particular, we study vacuum excitation probabilities 
(VEP), i.e., probabilities To,g->e of the detector evolving 
from its ground state \g) to its excited state je) when the 
held is initially in the vacuum state |0). We calculate 
the vacuum excitation probability or vacuum response for 
the four different detector models under consideration, 
up to leading order in perturbation theory. Calculation 
of the VEP serves a double purpose: (a) to demonstrate 
the need for renormalization of the models 2, 3 and 4 
in Table I and, after successful renormalization, (b) as a 
means to gauge whether these two models can be viewed 
as equivalent. 

For concrete calculations, we consider both bosonic 
and fermionic quantum helds F quantized in the cylinder 
with toroidal spatial sections consisting of the following 
product of intervals with periodic boundary conditions: 

B = [-L/2, V 2 ]" C R” . 


(17) 
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Therefore, 

F{t,x + 8) = F{t,x) V^eL-Z”. (18) 

The zero mode. — In the cylinder quantization, the field 
boundary conditions allow for constant classical solutions, 
that is, solutions of vanishing momentum k and energy 
Wfc. These degrees of freedom mathematically behave 
like free particles, rather than harmonic oscillators. They 
therefore have to be quantized differently from the rest of 
the field modes. Since there is no canonical vacuum state 
for such a zero mode, further subtleties appear when we 
couple the zero mode to particle detectors [40]. 

One could get rid of all the problems of the presence 
of a zero mode by using alternative boundary conditions 
which, unlike the periodic boundary conditions, do not 
permit solutions of zero momentum. A simple such choice 
is the use of Dirichlet boundary conditions, F{t, a;) = 0 for 
all X € dB, such that solutions of zero momentum vanish 
everywhere. These boundary conditions are straightfor¬ 
ward to implement for scalar fields, but full Dirichlet 
boundary conditions are too restrictive for spinor fields 
in the Dirac representation on C"* because the only solu¬ 
tion would be = 0. Instead, a weaker variant can be 
realized: imposing Dirichlet boundary conditions on the 
upper (or lower) two components of a four component 
spinor is possible; see for example [41]. However, it is not 
as clear how to construct a basis of the solution space that 
consists of eigenfunction of (corresponding to energy 
eigenfunctions in quantum theory) and <73 (corresponding 
to spin eigenfunctions), which makes quantization of the 
theory too involved. In this work we will therefore keep 
to periodic boundary conditions and simply exclude the 
zero mode from interactions, as described above. 

Fortunately, in most cases, it is possible to choose “safe” 
states for the zero mode that would minimize its impact 
on the dynamics of the particle detectors [40]. Based on 
this fact, in this work we will assume ad hoc that the 
detector does not couple to the zero mode of the field, as 
it is often done in quantum optics. The question in what 
circumstances dropping the zero mode can be justified 
is not vital to the discussion here, and is still subject to 
active research. 

Defining the vacuum response. — We assume that both 
field and detector are in their respective free ground states 
at time to. The state of the composite system is then 
characterized by the density operator 

pito) = \0,g;to) {0,g;to\ = \0,g) {0,g\ . (19) 

The coupled system is allowed to evolve unitarily until 
time t = to+T; time-evolution is performed using the time- 
evolution operator U generated by the total Hamiltonian 

H = Hf + Ho + idint , (20) 


and the probability for the detector to be excited is en¬ 
coded in the corresponding component of the density 
operator, after tracing out the field: 

do,3->e(I, to) = {e\trF pit)\e) . (22) 

At the end, we take the limit 

Po.g-S-e = lim Po,g^e{t,to) (23) 

to —^ —oo 
t —>-+oo 

such that the actual duration of the interaction is de¬ 
termined by the switching function x(t) alone: at times 
outside the support of y, both systems evolve freely. 

It is convenient to reformulate Eq. (22) in terms of 
state vectors in the interaction picture: 

Po,g^e{t,to) = '^\{a,e;t\U{t,to)\0,g-,to)f , (24) 

a 

where |a) is an orthonormal basis of the Hilbert space 
of states of the free quantum field F. Recall that the 
interaction picture is related to the Schrodinger picture 
through the partial inverse time evolution Uo{t) = 
which only takes into account the free Hamiltonian Ho = 
Hf + Ho [42]. 

Finally, we set up the usual perturbation theory under 
the assumption that A is small. To this end, the time- 
evolution operator is expanded in the Dyson series 

U{t,to) = J2 -^U^"'’it’to) , (25) 

n! 

n—O 

where the operator at order n is [43] 

1 ^ rt 

U‘^^\t,to) =—TH dUH,M ■ (26) 


For n = 0 we have = 1 (no change over time at 
all), but since {e\g) = 0, order zero does not contribute 
to do.g->e in Eq- (24). Therefore, 


do,(?—(^5 Iq) ~ -iA {a, e;t\U''^\t,to)\0, g;to) 

a 



(27) 


to leading order in A, where is simply the probability 
amplitude for the system to evolve from |0,(/) to ja, e) at 
leading order. 

In the following, we calculate the VEP for quantized 
scalar fields in linear coupling, as well as for quantized 
scalar and spinor fields in quadratic coupling (all models 
in Table I). 


A. Linear coupling to quantized scalar fields 


where Hf is the Hamilton operator of the free field. Ho 
is given by Eq. (5), and i7i„t depends on the model under 
consideration. The system ends up in the state 

p{t) = U{t,to)p (t, to) (21) 


1. Vacuum excitation probability 

In preparation for the derivation of new results con¬ 
cerning fermionic detector models, we begin by rederiving 
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the standard case of a quantized real field $ (see, among 
others, [34, 35, 44]). In order to work as generally as 
possible, we will avoid choosing a particular spatial pro¬ 
file and switching function. We couple an UDW-type 
detector to the quantized scalar field using the interaction 
Hamiltonian Eq. (14) (model 1 in Table I). We are dealing 
with a detector at rest, so its trajectory is 

x{t) = Xq , t(r) = T . (28) 

The time-evolution operator at leading order reads 




/ ‘^yp{xo,y)^{t\y) (29) 


’to 


/R" 


according to Eq. (26). Notice that for a sufficiently local¬ 
ized spatial profile p (as for example a Gaussian profile), 
if the localization length a is much smaller than the 
compactification length L, we can extend the integration 
region to the full space to a good approximation. 

Here, the monopole operator in the interaction picture 
is 


Solving the Klein-Gordon equation in a cavity with pe¬ 
riodic boundary conditions yields (see Appendix A) for 
the spatial part of the mode function 

p,,{x) = , (34) 

so the final expression for the VEP reads 



Z. Dependence on the dimension of spacetime 

The momentum sum in Eq. (35) depends on the spatial 
dimension n, since k = (fc^, ...,/c"). Therefore the 

convergence behavior of Po,g->e a-lso strongly depends 
on the dimension. To demonstrate this, we assume the 
detector to be pointlike. 


p(t) = (T- . (30) 


p{x,y) = S{x - y) , (36) 


Similarly, the field operator 




v^2a^ L 




flfce ifkix) 

^ ifl{x)\ , (31) 


where ujk is the energy of a single mode, the ladder opera¬ 
tors a/c satisfy the usual canonical commutation relations 
Eq. (A14), and the mode functions pk Eq. (AlO) form 
an orthonormal basis of the solution space of the Klein- 
Gordon equation. The field is restricted to a cylinder, so 
the momentum spectrum is discrete. 

Since the detector does not couple to the zero mode, 
we have 

O'jr 

fe G — • Z" \ {0} . (32) 

Lj 

Note that we have not yet chosen a particular form for 
the switching function or spatial profile. 

From Eq. (27), the only nonvanishing amplitudes are 
A[^\. Therefore the time evolved states will be superpo¬ 
sitions of the ground state |0,g) and states of the form 
|lfe, e), which feature a single quantum in the field. Sum¬ 
ming up the squared moduli of the amplitudes and taking 
the time limits leads to the VEP: 


P< 


0,g->e 




2 OJk 


p{xo,y) Pk{y)dy 


X{t) dt 


■ (33) 


and switch it on and off abruptly (sudden switching) by 
means of a window function: 


x{t) 


1 if t G [to, to + T] 
0 else 


(37) 


This is the original UDW detector introduced by DeWitt 
[3]. Plugging y and p{y,-) in Eq. (35) and performing 
the time integrals, the VEP simplifies to 


Pr 


2A 


O.g-^e 


Zi/\ V 

~ rn 


1 


UJk{Xl + UJk)' 


sm 


n 




T] . (38) 


In the massless case, the momentum sums are made ex¬ 
plicit by writing out the dispersion relation 

OJk = |fc| = + I 2 + • • • + In , (39) 

where G Z \ {0}. 

In one spatial dimension, we have 


P 


0,g^e 





h 



sin^ 


TT 

L 




This sum is bounded by 


S 3 2 ^ 3 

/i = l 1 


( 40 ) 
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which is convergent. 

For dimension n = 2, the probability is 



which is similarly bounded by the convergent double sum 


In higher dimensions convergence can only get worse, 
because even more sums will appear. For any dimension 
(l,n) with n > 3, the VEP is divergent [35]. 

The above results show that the VEP strongly depends 
on the dimension of the underlying spacetime. In the 
above example, the vacuum response of a pointlike detec¬ 
tor with sudden switching, coupled to a quantized real 
field is convergent in (1,1) and (1,2) dimensions, but 
divergent in all higher dimensions, necessitating regular¬ 
ization. 

3. Regularization through detector profile 
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In n = 3 dimensions, however, the VEP diverges. Its 
formal expression is 
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As Q,L{2tt) ^ > 0, we can estimate 
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By exploiting cos(x) = sin(x-I-7r/2) and the Pythagorean 
trigonometric identity, it is straightforward to prove that 
the right-hand side converges for all choices of parameters 
if and only if 


Those divergences of the VEP can be understood as a 
result of the pointlike structure of the detector as well as 
the sudden switching [35]. In the linear coupling case, the 
divergences can be regularized by adiabatically switching 
the detector, or by “smearing” the interaction between 
detector and field over the spatial profile of the detector. 
Let us study both possibilities in two separate examples. 

Gaussian switching and pointlike detector .— Consider 
the probability Eq. (35), still with a pointlike detector 
Eq. (36), but this time employing a Gaussian switching 
function. 


X{t) = exp 



(45) 


which is known to regularize the response of a resting 
detector in (1,3)-dimensional Minkowski spacetime [34]. 
We generalize this result to (l,n) dimensions: the time 
integral is readily solved, giving 


g—t^/(2T^) gi(n-i-wfc)tg|^ _ yfi2TrT‘^ g— 


Substituting this in Eq. (35) yields 


(46) 
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which is bounded by the convergent sum 
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does, where A = f7P(27r) This sum is in turn bounded 

by 


EEE [A + h+h+hr^ (44) 
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This shows that a suitably smooth switching function is 
able to regularize the leading order excitation probability 
in arbitrary dimensions. 

The excitation probability in this example is finite, but 
nonzero, even though the detector was in the ground state 
and the field in the vacuum. This should not be surprising 
considering that we have a time-dependent Hamiltonian: 
When the interaction is switched on, the Hamiltonian of 
the composite system is changed, and the state jO,;/) is 
no longer an energy eigenstate. In consequence, the state 
begins to evolve nontrivially, and there is a finite proba¬ 
bility for the detector to be measured in its excited state 


which diverges logarithmically. 



|e). The more adiabatically the interaction is switched 
on, the lower is the probability to excite the detector [34]. 
This is reflected in Eq. (47): the larger T, that is, the 
slower the detector is switched, the smaller Pg^g^e- In 
the limit T —>■ oo, corresponding to adiabatic switching, 
Po,g^e —> 0 as expected. 

Sudden switching and Gaussian detector profile .— Now 
consider a Gaussian detector profile. 


p{xo,y) = i2Tra‘^) exp ^ 

together with a sudden switching. In this case: 
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in the literature, namely, model 4 in Table I pioneered by 
Takagi in 1985 [33] and extensively used, e.g., in [20, 36" 
39]. We will assess the physicality of this class of models 
by considering the vacuum response of spatially smeared 
detectors. 

Let be a spinor field, that is, a field taking values in 
some representation space C™ of the complexified Clifford 
algebra Cliff (l,n)®. In order to be able to do explicit 
calculations, we choose the irreducible Dirac representa¬ 
tion of Cliff(I,3) on in 4 spacetime dimensions. In 2 
spacetime dimensions, a related but reducible represen¬ 
tation of Cliff(l, I) on C"* is used, which yields similar 
spinor mode functions. Details on these conventions and 
standard results in classical and quantum field theory can 
be found in Appendix A. 

We will work again in the interaction picture, where 
the operator of the quantized field is 

'i>{t,x) = 

+ &I,se‘“'“*V'fc,s-(®) • (54) 

/c,s 


and therefore again convergent in arbitrary dimensions— 
the spatial smearing of the interaction between detector 
and field is also able to regularize the vacuum response 
of the detector. Moreover, Po.g-ye 0, if the detector is 
completely delocalized, that is for u —> oo. 

General effect of spacetime profile .— The two preceding 
examples show that it is possible to regularize the VEP to 
leading order in perturbation theory. But which combina¬ 
tions of switching function y and spatial smearing p are 
able to regularize the VEP? To address this question, it is 
helpful to use the detector’s spacetime profile /: the VEP 
of the quantized real field in Eq. (35) can be reformulated 
as 


The ladder operators ak,s, bk,s satisfy the canonical anti¬ 
commutation relations Eq. (A60). The spinor-valued 
mode functions tpk.s.e span a solution space of the 
Dirac equation. Depending on dimension and mass, the 
mode functions take slightly different forms, namely, (a) 
Eqs. (A28) and (A37) for massive and massless fields in 
(1,3) dimensions; and (b) Eqs. (A28) and (A37) in (1,1) 
dimensions. 

This quantized spinor field is coupled quadratically to 
a resting UDW-type detector via the interaction Hamilto¬ 
nian Eq. (15), which is model 4 in Table I. Using Eq. (26) 
we then obtain the leading order contribution to the time 
evolution operator: 


P 


O.g^e 
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^ — /(Wfc -k D, k) 
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(52) 


where / is the (1, n)-dimensional Fourier transform of the 
spacetime profile /: 

/(w$,|)=/ dx [ df/(a;,t)e"’“«‘e'^"’ . (53) 

J J — oo 

Therefore, a given spacetime profile does regularize the 
VEP to leading order in (l,n) dimensions if and only 
if the modulus squared of its (1, n)-dimensional Fourier 
transform / decays fast enough in the UV such that the 
n-fold sum in the above expression is finite. 


B. Quadratic coupling to quantized spinor fields 

We now turn to the pivotal question: How can we 
construct a similar particle detector model for quantized 
spinor fields? To this end, let us begin by investigating 
and extending a model that has been commonly employed 


U^^\t,tfi)= [ dt'xit') p{t') 
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The vacuum response Eq. (27) yields 
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There are obvious differences to the vacuum response 
of the usual UDW model (model 1 in Table I, i.e., a 
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monopole detector linearly coupled to a quantized scalar 
field). Comparing Eq. (33) with Eq. (55), at leading order 
in A, there are two kinds of terms instead of one. 

The first term corresponds to excitation of the detector 
by creation and subsequent annihilation of a field quan¬ 
tum, with amplitude More precisely, the detector is 

excited by creating and annihilating an antiparticle from 
the vacuum. However, the equivalent process featuring a 
particle does not occur. 

The second term represents processes where the detec¬ 
tor is excited by emission of a particle and an antiparticle, 
with amplitude t • These stood to be expected 

in the light of our discussion regarding fermion number 
conservation in Sec. IIB. Note that the momenta of the 
particle and antiparticle are not related since the detector 
is “heavy” and can absorb any amount of momentum. 

We can make the mode functions ipk,s,e explicit using 
Eqs. (A28), (A37), (A45) and (A53). In all four cases 
considered here, massive and massless helds in (1,3) and 
(1,1) dimensions, the VEP when coupling quadratically 
to quantized spinor fields can be brought into the form 


Pf 


0,g->e 


+ E 


A2 




2L2 

(wfc -I- m){ujp -I- m) 




k 


P 


ujk+Tn ojp + m 


p{xo,y)e *(''+P)3^dy 


/R" 



r+oo 

2 ■ 

X 

/ x(t) +“?>)* dt 



J —OO 



(57) 


where n = 1 or n = 3. The case of a massless field 
is correctly recovered by setting m = 0, such that the 
first term vanishes. Notice that there is an additional 
problem with the first term in Eq. (57) for m ^ 0. Namely, 
the term is proportional to Even in the least 

divergent scenario, i.e. n = 1, this sum diverges like 

This divergency is fundamentally worse behaved than 
those encountered in the standard bosonic Unruh-DeWitt 
model, which arise from the pointlike structure of the 
detector [13, 20], or its overly abrupt switching [34, 35]. 
The first term in Eq. (57) cannot be regularized by way of 
a smooth spacetime profile, and therefore the divergency 
in model 4 cannot be cured via this common method. This 
is because the spatial integral containing the (normalized) 
detector profile factors out and yields a factor of 1. The 
time integral here is not momentum dependent and merely 
contributes an overall factor which cannot regularize the 
divergent sum. 

Comparing with the Unruh-DeWitt model, there are 
three conceivable origins of this new divergency: (a) it 
could be algebraic on the level of ladder operators, because 
we are now dealing with a field obeying Fermi statistics 
instead of Bose statistics, (b) it could be analytic due to 


the field’s spinorial structure and different inner product 
compared to the scalar fields, (c) or the origin could be 
the coupling that is now quadratic in the field instead of 
linear. The question is easily settled by computing the 
VEP of a quantized complex field coupled quadratically 
to an UDW-type detector, keeping properties (a) and (b) 
the same as in the usual Unruh-DeWitt model. 


C. Quadratic coupling to quantized scalar fields 


Let $ be a quantized real or complex field (models 2 
and 3 in Table I, respectively), coupled quadratically to 
a resting UDW-type detector via the interaction Hamilto¬ 
nian Eq. (16). Inserting the interaction Hamiltonian in 
Eq. (26) yields 


U^^\tM) = 

[ dt'x{t')p{t') [ dyp{xo,y)d>\t',y)<P{t',y) , (58) 
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where the field operator (in the interaction picture) is 
given by Eq. (31) in case of a quantized real field. If the 
field is complex 


+ pUx)] ■ (59) 


Calculation of the VEP Eq. (27) yields 
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(60) 


in both cases at leading order. Notice that this is no longer 
true at higher orders. This expression has two terms like 
Eq. (56), representing analogue processes: The first term 
again corresponds to excitation of the detector by emission 
and reabsorption of an antiparticle with amplitude A^l, 
and the second term corresponds to processes where the 
detector is excited by emission of a particle-antiparticle 
pair with amplitude A[^\ . Explicitly inserting the 

mode functions Eq. (34) for periodic boundary conditions 
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gives 
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(61) 


The first term is proportional to thus diver¬ 

gent in the same way as for spinor fields in Eq. (57). Since 
this divergence appears in exactly the same form for both 
scalar and spinor fields, its origin is clearly the quadratic 
coupling. Also note that for quantized scalar fields the 
divergence arises even if to = 0. 

It is not surprising that these divergencies arise: In¬ 
teractions containing second or higher powers of fields 
tend to lead to divergencies in observables, which require 
renormalization. In particular, the spinor field detector 
model by lyler and Kumar [32] faces the same problem. 


IV. RENORMALIZATION OF QUADRATIC 
MODELS 

Before investigating the renormalization of the quadrat- 
ically coupled detector models Eqs. (15) and (16), let us 
briefly review the literature. 


A. Literature review 

In Ref. [17], Hinton introduces an UDW-type particle 
detector model coupled to a quantized real field $ through 
the interaction Hamiltonian 

i7i„t = , (62) 

which was the template for Eq. (16). Hinton calculates 
the VEP in [17, eq. (14)[, remarks that the expression is 
formally divergent and requires renormalization but no 
explicit renormalization scheme is given. 

Takagi, who suggested coupling a quantized spinor 
field quadratically to DeWitt’s two-level detector through 
Eq. (2), followed a different approach. In the original 
publications, see Refs. [20, 33], Takagi points out an in¬ 
finite term in the VEP at leading order but argues that 
this term can be dropped. Later publications employing 
Takagi’s detector model or drawing on his results implic¬ 
itly follow the same reasoning [36-39]. A key point is 
that instead of the total excitation probability, excitation 
rates are considered. As is well known, the calculations 
of rates involves fewer integrations and are therefore gen¬ 
erally more regular. In addition, Takagi regularizes the 
detector response by using a smooth, namely exponential 


switching function. As a consequence, in the excitation 
rate 

R= lim Po,g^e{t,to) 

to—> —oo t — to 

t —>-+oo 

the divergence in Po,g->e is offset by a divergence in the 
denominator. 

Takagi’s expression for the rate [20, eqs. (8.5.2) to 
(8.5.4)] can easily be obtained in the notation developed 
earlier by using Eqs. (15), (26) and (27) in Eq. (63) and 
by using a pointlike detector with exponential switching 
function: 

p{x, y)=5{x-y) , xi.t) = (64) 

The limit s —>■ 0 is taken at the end of the calculation in 
Ref. [20] so that ultimately the detector is switched on 
and off infinitely slowly. One obtains the rate 

R = \^\{e\p\g)\^F{n) . (65) 

Here, 


E(Q) = lim lim [ dti [ dt; 

s—^Oto^ — oo 1. 1. 


‘t-to 


X ^ ( 66 ) 


which is usually called the detector response function [14], 
and 


S{ti - t2) = (0|4'(ti,x)T(ti,x)T(t2,®)^'(t2,®)|0) (67) 


is a 4-point correlation function. This step is not made 
explicit in Ref. [20]; rather, Takagi directly argues that 
the exact form of the switching function is not relevant, 
since it has only been introduced as a way to regularize 
the detector response. Therefore, he argues, one is allowed 
to temporarily replace the exponential switching function 
with a Gaussian one, and to exploit that only the limits 
t —>■ oo, to —oo and s —)■ 0 are of interest—see end of 
section 3.3 in Ref. [20]. Finally, the Gaussian is replaced 
ad hoc with an exponential function again. If we accept 
Takagi’s reasoning, the above response function can be 
rewritten as 

/ OO 

dte"‘^(‘)e""l‘IS'(f) , (68) 


which indeed corresponds to [20, eqs. (8.5.2) to (8.5.4)] 
(with two differences: in the original paper, the factor 
due to the switching function is suppressed, and 
A = 1). In this notation, the divergence we encountered 
earlier is now hidden in the correlator S: by Wick’s 
theorem it can be expanded as 

S{ti-t2) = [trS'“(0)]^-|-tr[S'“(ti-t2)5''''(G-f2)] (69) 
where S~, S'^ are the Wightman functions 

S~{ti,t2)i = (0|TB(t2,a;)«'^(fi,a;)|0) 
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and the first term is infinite, as expected: 

tr 5'“(0) —)■ oo . (71) 

Takagi treats this infinite term as formally constant in 
Eq. (69), which allows one to obtain a delta distribution 
from the time integral: 


Fin) = Sin) [^^-(O)]^ + lim / dt'e-‘^(*')e-"l‘'l 

X tr S~it 2 ,ti)S^iti,t 2 ) ■ (72) 

He concludes, therefore, that the term trS'“(0) can be 
dropped, as long as the detector gap is finite, O > 0, 
thereby obtaining a finite excitation rate R. 

This argumentation is not entirely satisfactory. In par¬ 
ticular, it merely yields a finite excitation rate. This is 
not enough to fully characterize the response of a particle 
detector. For that, one would need to be able to com¬ 
pute its full density matrix, which requires computing 
probability amplitudes, not rates. Takagi’s techniques, 
however, do not avoid divergences in probabilities. The 
challenge remains to fully renormalize particle detector 
models with an interaction Hamiltonian that is quadratic 
in the field. 


B. Renormalization at leading order 

As we discussed above, the first term in Eqs. (57) 
and (61) contains a divergence that cannot be regular¬ 
ized by means of a spacetime profile, unlike in the usual 
linearly coupled UDW model. As we will show now, the 
occurrence of such divergencies is not surprising from the 
point of view of a quantum field theoretical renormaliza¬ 
tion theory, and can be handled straightforwardly. To this 
end, let us compare the two problematic detector mod¬ 
els to quantum electrodynamics (QED). The interaction 
Hamiltonian of QED, Hqed oc is also quadratic 

in the spinor field (electron field), the only difference being 
that the coupling is not to a scalar detector but to the vec¬ 
torial photon field. Indeed, QED has a similar divergence, 
namely the well-known divergent tadppole subdiagram. 
In QED, it is renormalized straightforardly by normal¬ 
ordering the interaction Hamiltonian (see, e.g. [45]). In 
the following two sections we demonstrate that and how 
the divergencies in the detector models are related to tad¬ 
pole diagrams, and we show that normal-ordering of the 
interaction Hamiltonian renormalizes the detector models 
as well. 


1. Structure of the divergence 

As mentioned before, there are two types of perturba¬ 
tive processes that excite the detector at leading order, 
thereby contributing to the VEP for quantized spinor 
fields in Eq. (56), or quantized scalar fields in Eq. (60). 



(a) Divergent first term. 



(b) Second term. 


Figure 1. Schematic representation of the processes leading 
to the state Eq. (73). The wiggly line stands for the detec¬ 
tor (dashed when in ground state, solid if excited), the solid 
straight line for the quantized scalar or spinor field. The time 
direction is from left to right. 


To see this, recall that at leading order, time evolution 
amounts to applying the interaction Hamiltonian to the 
initial state |0, g) and integrating over the entire duration 
of the interaction. In case of a spinor field, the final state 
is Given the quadratic coupling in the interaction 

Hamiltonian ~ Tlk for the fermionic case and ~ <I>1<I> for 
the bosonic case, the algebraic structure of the leading 
order time evolved field state is of the form: 

IV') ,!<(') = E |0) + E ip) • (^3) 

k k,p 

The spinorial degrees of freedom are suppressed and Ek, 
Fif: p are coefficients depending on the momenta (and for 
spinor fields additionally on spins). In other words: When 
interacting with the field in the vacuum state, the detector 
can be excited (to leading order in perturbation theory) 
either by emission and subsequent annihilation of an an¬ 
tiparticle [divergent first term in Eqs. (56), (60) and (73)], 
or by emitting a particle-antiparticle pair (second term in 
the same equations). Since we are only interested in the 
state of the detector after the interaction, both types of 
processes contribute to the vacuum excitation probability. 
They are visualized as diagrams in Fig. 1. Note that these 
simplified diagrams (used for illustration) are not strictly 
Feynman diagrams since the detector is not yet second 
quantized, as we will see in more detail in Sec. VHI when 
we build the detector-field interaction Feynman rules. 

The following discussion is limited to the case of a 
quantized spinor field, but runs analogously for quantized 
scalar fields. In comparing the spinor field detector with 
QED, the spinor field corresponds to the electron field 
in QED, and the detector to the photon field. Thus, the 
detector ground state translates to the absence of photons, 
while the excited state corresponds to the presence of one 
photon. In this sense, the processes in Fig. 1 are analogous 
to the QED diagrams drawn in Fig. 2. In particular, the 
problematic first term in Fig. 1 corresponds to the tadpole 
diagram in Fig. 2. The divergence of its amplitude is 
mirrored by QED: the amplitude of the tadpole diagram 
is divergent as well; see ]45]. 

While the “heavy” detector guarantees energy and mo¬ 
mentum conservation in Fig. 1, both QED processes in 
Fig. 2 are dynamically forbidden. They may, however, 
appear as subdiagrams in larger processes, and it is possi- 
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(a) Tadpole diagram (b) QED vertex 


Figure 2. QED Feynman diagrams corresponding to the pro¬ 
cesses in Fig. 1. 


ble to calculate their contribution to the total amplitudes. 
Indeed, the analogy between the left-hand diagrams in 
Figs. 1 and 2 extends to their respective amplitudes: ac¬ 
cording to Eqs. (27) and (55), the amplitude of the left 
diagram in Fig. 1 is 

4^1 = -iAe-‘^* 

x(0,e| [ x(^')M(^')^(^'.ai)^'(t',ai)dt'| 0 , 5 ) 

Jto 

^ (74) 

when suppressing the time integration. The detector was 
assumed to be pointlike for simplicity, since we know 
already from Eq. (56) that the spatial profile has no 
regularizing effect. Using the fermionic time-ordering 
operator, this can be rewritten as 

(0| TA(t,£c) a;)|0) 

= -(0|TT^(t,ai)^.4(t,ai)|0) , (75) 

which is essentially the Feynman propagator: 

Sf{x - y)i = (0|T4/^(a:)TB(j/)|0) . (76) 

We find 

~ trS'i?(0) —>■ oo , (77) 

where the trace runs over the spinor indices. 

By applying position space Feynman rules to the tad¬ 
pole diagram in Fig. 2, on the other hand, one obtains 
the amplitude 

tr[7^S';’(0)] —>■ oo , (78) 

(see e.g. [45]). The divergence of the two amplitudes 
Eqs. (77) and (78) has a common structure: they di¬ 
verge because the Feynman propagator of the quantized 
spinor field is ill defined when closed in itself. The gamma 
matrices in the second amplitude merely appear because 
electrodynamics has a vector coupling Eq. (3), as opposed 
to the scalar coupling Eq. (15) for the detector model. 

In conclusion, it is justified to interpret the divergence 
found in the quadratically coupling detector models as 
the analogue of tadpole diagrams known from quantum 
field theories such as QED. 


2. Renormalization by normal-ordering 

Tadpoles are renormalized by normal-ordering of the 
interaction Hamiltonian. This amounts to setting the am¬ 
plitude of the tadpole to zero, or—equivalently—ignoring 
any diagram containing tadpoles as subdiagrams (see, for 
instance, [45]). The same procedure works for both the 
detector models as well: normal-ordering the interaction 
Hamiltonians Eq. (15) leads to 

AlW = -iAe-‘^* 

X (0, e| [ x(i') R(i') ■ x) 4'(t', x): dt' |0, g) 

Jto 

= 0 (79) 

instead of Eqs. (74) and (77), such that the divergence is 
renormalized to zero. Equivalently, for detector models 
coupled quadratically to quantized real or complex helds 
via the Hamiltonian Eq. (16), 

AlW~(0|:<I>t<i>:|0) = 0. (80) 

The above renormalization can be interpreted as shift¬ 
ing the expectation value of the interaction energy from 
infinity to zero: If detector and field do not interact, |0,g) 
is the ground state, with energy eigenvalue zero: 

(0,5|i7|0,5)=0. (81) 

Here, H = for quantized scalar fields, and 

H = Hd + H<s, for quantized spinor fields as given in 
Eqs. (5), (A15), (A19) and (A61). However, as soon as 
the coupling is switched on, A > 0, we would naively have 

(0,5|i7 +Flint |0, g) ^oo (82) 

for the interactions quadratic in the field Eqs. (15) 
and (16). This would be unphysical and normal-ordering 
the interaction Hamiltonian is therefore appropriate to 
ensure that 

{0,g\H + :Hi^t--\0,9) = 0 ■ (83) 

In comparison, for linear coupling Eq. (14) the energy 
expectation value vanishes automatically. 

C. The renormalized models and their VEP 

In summary, renormalization of the spinor field detector, 
model 4 in Table I, is achieved by coupling through the 
Hamiltonian 

= Xxit) g f p{x{t), y) : T(y)5'(y): dy . (84) 

JR" 

Detector models 2 and 3 for real or complex fields are 
similarly renormalized by coupling through 

= \x{t)y [ p{x{t),y) :$'l'(y)4>(y): dy . (85) 
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The vacuum response of the renormalized spinor field 
detector is obtained by dropping the first term in Eq. (56): 


Pf 
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Explicitly, for periodic boundary conditions: 
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Similarly, the leading-order vacuum response of the renor¬ 
malized real and complex field detector becomes 
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(89) 


by dropping the first term from Eqs. (60) and (61). 

As will be demonstrated in the next section, the VEPs 
Eqs. (87) and (89) may still be divergent even after this 
renormalization and further regularization of the VEP 
is necessary. Recall that the VEP of detector model 1 
coupling to a quantized scalar field could be rewritten 
in terms of the Eourier transform of its spacetime profile 
/ in Eq. (52). The same is possible for the quadrati- 
cally coupling detectors: using the Fourier transformation 
Eq. (53), the VEP Eq. (89) for quantized scalar fields can 
be rewritten as 


The probability in the case of Eq. (87) is 


P 


0,g->e 


A^ {ojk + m){up + m) 

2T^” WfcWp 


X 


Wfc 


p 

ujp + m 


f{-d-k-p) 


(91) 


Thus, the UV behavior of its Fourier transform again 
decides which spacetime profile / is suitable for regular¬ 
ization. Note that the prefactor in front of the / decays 
faster in the UV for quantized scalar fields than for quan¬ 
tized spinor fields. This means that spinor field detectors 
require a “stronger” regularization. 


V. COMPARISON OF PARTICLE DETECTOR 
MODELS 

So far, we have discussed how to obtain a finite vacuum 
response through renormalization and regularization for 
four different detector models. We will now investigate 
whether any of the three scalar field detector models 
(1,2 and 3 in Table I) is comparable to the spinor field 
detector, in the sense that a comparison of their respective 
responses offers insight into the properties of the field they 
probe. 

For the different models considered, the relevant fea¬ 
tures that can have a significant impact on the response 
of UDW-type detectors for quantized scalar and spinor 
fields are 

1. The nature of the coupling: linear versus quadratic. 

2. Internal degrees of freedom of the field: for example, 
the C/(l) charge. 

3. The field statistics: bosonic or fermionic. 

4. The analytic structure of the field: scalar or spino- 
rial. 

From the results obtained in the previous sections, we 
discuss below each of these points individually. 


A. Linear vs quadratic coupling 

The coupling has a profound influence on the response 
of an UDW-type detector: The vacuum response of model 
1—i.e. coupling linearly to a quantized real field—given 
in Eq. (33), is fundamentally different from the response 
[Eq. (88)] of model 2 to the same field when using the 
(renormalized) quadratic coupling of Eq. (85). In order to 
assess the differences between the response of detectors 
to fermionic (spinor) fields and bosonic (scalar) fields, we 
will now compare models 2 to 4 that couple to the same 
power of the fields. 


B. Charged vs uncharged field 


(90) Quadratic coupling to quantized scalar fields gives the 

same VEP Eq. (88) for uncharged (i.e. real) and charged 
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(i.e. complex) quantized scalar fields to leading order in 
perturbation theory. However, as we shall see in Sec. VIII, 
the detector response is not identical in general; the field 
charge does have an influence. Since spinor fields are 
generally charged, one should always use a charged scalar 
field when comparing quantized scalar fields with quan¬ 
tized spinor fields in order to single out effects coming 
exclusively from the field statistics or analytic structure. 
In other words, only model 3 and model 4 can be rightfully 
compared. 


C. Boson vs fermion statistics 


The statistics of the probed field influence the reaction 
of an UDW-type detector in two ways. First of all, the 
Pauli exclusion principle will prevent the creation of more 
than one quantum of a given charge, momentum and 
spin in a fermionic field. This restricts the set of possible 
processes that, for example, excite the detector. This 
would not be the case for bosonic fields. Secondly, the 
anticommutation of fermionic operators implies that two 
distinct processes leading to the same final state may have 
a relative minus sign in their amplitudes. Fermionic fields, 
unlike bosonic fields, can therefore have cancellations 
between amplitudes. 

Nevertheless, the field statistics do not influence the 
VEP of the four detector models discussed here at leading 
order: merely a single quantum (linear coupling), or at 
most a particle-antiparticle pair (quadratic coupling) is 
created from the vacuum of the field, so the Pauli exclusion 
principle has no effect. And since there is only one possible 
process contributing to the VEP, no cancellation between 
different processes occurs. This implies that the leading 
order VEP remains unaffected by the field statistics. 

As soon as the field is initially not in its vacuum state, 
however, the Pauli exclusion principle does affect the 
response of a particle detector to a fermionic field. And 
even for the VEP, the exclusion principle and different sign 
amplitude interference will be relevant at higher orders in 
perturbation theory. Interestingly, in dimensions larger 
than two, the spin-statistics theorem ties together field 
statistics and the spin of a field (which is reflected in 
its analytic structure). This makes it generally hard to 
distinguish the effect of the analytic structure and the 
effect of field statistics. In the leading order vacuum 
response, however, due to the uniqueness of the leading 
order VEP described above, it is possible to independently 
study the effect of the spinor structure versus a scalar 
structure. 


D. Scalar vs spinor field 


Comparing Eqs. (86) and (88), one finds that the VEP 
of model 3 has summands proportional to 


4wfcWp 


p{xo,y) p*k{y) (p*{y) dy 


/R" 


p{xo, y) y) y) dy 


'R" 


, (92) 


while the VEP of the spinor field detector, model 4, has 
summands proportional to 



p{xo, y) '4’k,s,+ {y) '4’p,T-{y) dy 


2 



p{xq, y) ^k,s,+ (t, y) 7° y) dy 


2 

(93) 


The sum over all spins in the second expression is, of 
course, not present in the first one for quantized scalar 
fields. The factor (dwfeWp)”^ in the first expression is 
simply part of the normalization factor of the scalar mode 
functions Eq. (A5). This can be made explicit by rewrit¬ 
ing both of the above expressions in terms of the full, 
time-dependent mode functions Eqs. (A5) and (A25) (at 
arbitrary time t). In this form, it is obvious that the only 
difference between the VEP for scalar and spinor fields is 
due to the different mode functions. 

So does it make any difference at all whether the de¬ 
tector is coupled to a quantized scalar or a spinor field? 
The answer is yes: Comparing the VEPs Eq. (90) and 
Eq. (91) at the end of the previous chapter, we found 
that quantized scalar fields are better behaved for large 
frequencies since the UV behavior of the mode functions 
is different; the scalar mode functions generally have a 
factor ~ (wfc)”^/^ which the spinor mode functions do 
not: As a simple example, compare the mode functions 
of massless fields in (1,1) dimensions: For the quantized 
(real or complex) scalar field, the particle and antiparticle 
mode functions are 

[see Appendix A, Eqs. (AlO), (A12) and (A17)], while 
the corresponding spinor mode functions for particles and 
antiparticles (in this order) are 


'ipk,s-{t,x) 


1 


1 


is 

Sgn{k) 

sgn(A:) 

is 


g-|-i|fe|t g-ifca: 


(95) 


[see Appendix A, Eqs. (A25), (A45), (A53) and (A58)]. 
The spinorial parts do no scale with k, while the scalar 
field has an additional factor 
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This difference has notable implications: it causes the 
quantized scalar fields to exhibit better convergence prop¬ 
erties when summing over all momenta—the VEP for 
detectors coupling quadratically to quantized scalar fields— 
Eq. (90)—will in general converge better than the VEP 
for detectors coupling quadratically to quantized spinor 
fields—Eq. (91). This is a dramatic difference. Eor in¬ 
stance, there are switching functions and spatial profiles 
that give a finite VEP for a detector quadratically coupled 
to a quantized scalar field but yield a divergent VEP for 
quantized spinor fields. In Sec. VE, examples where this 
situation actually occurs are discussed. 

Superficially, the different normalization conditions, 
stemming from a mathematically different inner prod¬ 
uct, are the reason for this difference: The scalar mode 
functions are orthonormal with respect to 


massless fields to 


P( 


0,s->e 


_ A2 1 

4L2 \kp\ 

k,p^0 ' 


pixo,y)e '('=+p)«d|/ 


r-l-oo 


xit) e'(^+l'=l+IPl)‘dt 


• (98) 


On the other hand, for the spinor field detector (model 4 
in Table I), one similarly obtains 


Pr 


2A 

“ r2 


0,g->-e 


T2 


k,p>0 


p{xo,y)e 


(99) 


p{xo,y)e 


f-l-oo 


x{t) dt 


i'pk ) ‘Pp) — 


‘Pkix)[do(fip{x)] - [dQipk{x)](p*p{x) dx 

= Sk.p , (96) 


directly from Eq. (86) by plugging in the mode functions 
and spinors Eqs. (A45) and (A53). 


1. Sudden switching and pointlike detector 


while the spinor mode functions are normalized according 
to 


{‘4^kyS,e j '4^p,r^5) — 


'^k s e(^) V^Pj^a(^) — ^k,p ^e,S • (91^) 


The additional time derivative in the product for scalar 
functions generates a factor cok- Thus, the normalization 
factor of the scalar fields contains an additional factor 
(wfc)“^/^ compared to the normalization of the spinor 
fields. 

Ultimately, the reason is that the equation of motion 
of scalar fields, the Klein-Gordon equation, is of second 
order and features a double time derivative, while the 
equation of motion of the spinor field, the Dirac equation, 
is of first order, containing only a single time derivative 
(see Appendix B). 

In order to demonstrate the infiuence this has on the 
response of the different detector models, we compare 
their VEPs in three simple examples on (1,1) dimen¬ 
sional Minkowski spacetime for massless fields: (1) sudden 
switching with a pointlike detector, (2) Gaussian switch¬ 
ing with a pointlike detector, and (3) sudden switching 
with a Gaussian detector profile. 


E. Examples in (1,1) dimensions 


In the most basic scenario, the detector is pointlike 
[delta spatial profile as given in Eq. (36)], and it is sud¬ 
denly coupled to the field as given in Eq. (37). This is the 
same situation that was considered in Sec. Ill A 2 for the 
case of a quantized real field with linear coupling (model 

I). 

Quantized scalar field .— Plugging Eqs. (36) and (37) 
in Eq. (98), the resulting VEP for the quantized scalar 
field is 


^ i: i: 


it^uf^ihhi^ + ii + hy 


X sm 


TT / QL , , , 


( 100 ) 


This double sum is convergent: we can estimate 


47r^ 


Po,g^e ^ T T 


1 


■ 


^ 1^2 (ll + hY 

The sums on the right-hand side converge if and only if 


oo oo ^ 

^=EE—^ 


+ 12 Y 


does, where the summation starts at 2 instead of 1. We 
can now exploit that I 1 I 2 > h + h for all k > 2 to find 
the upper bound 


Let us start by recalling the general expressions for the 
massless (1,1) dimension scenario. Eor model 3 in Ta¬ 
ble I (quadratic coupling quantized scalar field detector), 
the VEP Eq. (89) simplifies in (1,1) dimensions and for 


^^EE 


1=212=2 


{ll + 12 ) 


where the sum on the right-hand side is convergent. 
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Quantized spinor field .— The formal expression for 
the VEP in case of a quantized spinor field is obtained 
by plugging Eqs. (36) and (37) in Eq. (99): 


, 2 °° °° 

^ E E 




X sm 


TT / QL , , , 


( 101 ) 


Similar to Eq. (43), it is easy to prove that this sum 
converges if and only if 


OO OO 


EE 


1 



2 


does. However, this sum diverges, since the corresponding 
integral is ill defined: 




dh 


1 

{A-\- li -\- l2f 


lim ln(a) — ln(2 + H) —>■ oo . 

a—¥oo 


Thus, the VEP Eq. (101) is divergent. 

This result nicely illustrates that the detector response 
to quantized scalar fields is better behaved in the UV: 
On one hand, all the scalar fields detectors (models 1- 
3) have a finite VEP in (1,1) dimensions, even if the 
detector is simply pointlike and uses sudden switching 
[see Eqs. (40) and (100)]. On the other hand, in the case 
of a quantized spinor field the VEP is already divergent 
in (1,1) dimensions for this spacetime profile. 


2. Gaussian switching and pointlike detector 


It was demonstrated in Sec. Ill A 3 that the original 
UDW detector (model I), which has a divergent VEP 
in dimensions (I,n) for n > 3, can be regularized by 
introducing a Gaussian switching function Eq. (45) of 
width T. We can prove that Gaussian switching also 
regularizes the spinor field detector, model 4, in (1,1) 
dimensions. 

Quantized spinor field .— Inserting Eq. (45) in Eq. (99), 
the VEP for a spinor field reads 


Pr 


0 , 3 ->e 


L2 


EEe-"^( q§-+h+i2Y ^ (^02) 

ll — l ^2 — 1 


This excitation probability is now clearly finite: the sum 
converges if and only if 


Quantized scalar field .— For model 3, the same switch¬ 
ing function yields 


2 \2 


Pr 
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T^X 

2tt 


^1^2 

lo = l 




(103) 


which is of course also finite since the summand decays 
faster than in Eq. (102) [and the VEP had been convergent 
in Eq. (100) even without regularization]. 

Note that for both field types, Po,g_>.e —> 0 in the 
adiabatic limit P —> oo, iust like for the linear coupling 
in Eq. (47). 


3. Sudden switching and Gaussian detector profile 


The third example discussed in Sec. Ill A 3 for model 1 
demonstrated regularization through a Gaussian detector 
profile Eq. (49) with variance cr > 0 and sudden switching. 
We repeat the example in the present case. 

Quantized spinor field .— In case of the quantized spinor 
field, the probability is 


Pf 


4A2 “ “ 


O.g-^e 
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X sm 
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(104) 


and convergence is harder to assess than for the Gaussian 
switching function in Eq. (102): the sine squared is posi¬ 
tive, bounded from above by 1 and can thus be ignored. 
For li = I 2 , the exponential factor is always 1, but then 
the denominator decays as as it happens, this is just 
fast enough. As a first step, we can estimate 


4A2 


Pf 
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OO OO 
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where A = ^ > 0. This sum converges if the 

corresponding integral 

g —^(^1—^ 2 )^ 

{h + 12 ) 



does. By substituting I = I 2 — h, 


I = 
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-Al 
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(1 + 2 ^ 1 )^ 




d;2 


does, where A = Att^T'^L ^ and B = VlL{2'it) and the 
latter expression is bounded from above by 



and since for I > 1 — li we have (^ -I- 2/i) ^ < l^^ 
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The VEP Eq. (104) is finite. 

Quantized scalar field. — For the quantized scalar field, 
the probability 
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is once again convergent because Eq. (100) already was 
even without regularization. 

Again, Po,g->e —> 0 for both fields if the detector is 
maximally delocalized by cr —)■ 00 . 


F. Comparability of the models 


From the analysis above we conclude that the differ¬ 
ences in the response of models 3 and 4 are only due to 
the analytic structure of the fields and the field statistics. 

The field statistics do not enter the leading order VEP, 
but will have a profound influence on the detector response 
at higher orders, or also at leading order if the field is 
initially not in the vacuum state. On the level of the 
analytic structure, model 4 is sensitive to the additional 
spin degree of freedom, and, more importantly, requires 
stronger regularization since it displays a worse-behaved 
UV response. Interestingly, it is not enough to regularize 
with a smooth switching function as opposed to the scalar 
case. Rather, the spinor model requires a spatial smearing 
in order to yield finite VEP. 

Note that only a qualitative comparison is feasible. A 
naive quantitative comparison fails because the coupling 
strength A has different dimensions in model 3 and 4: 
Scalar fields on (1, n)-dimensional spacetime have mass 
dimension (n — l)/2, while spinor fields have mass dimen¬ 
sion n/2. To obtain a Hamiltonian of mass dimension 
1, the respective coupling constants therefore need to be 
of different dimensions. A direct comparison is therefore 
only possible with regards to the functional dependence 
on model parameters like the detector trajectory x(t), the 
detector gap H, and more concretely the interaction time 
T and the detector size a. In this sense, the leading-order 
vacuum response of models 3 and 4 is equivalent in all 
(convergent) examples discussed in Sec. VE. 

In conclusion, for finite size detectors after renormaliza¬ 
tion, the pair of models, model 3 (for quantized complex 
fields) and model 4 (for quantized spinor fields) can be 
used to reliably study the differences between fermionic 
and bosonic fields via particle detectors, at least at leading 
order in perturbation theory. 


VI. COMPUTATION METHODS AT 
ARBITRARY ORDER 

Up to now, the discussion was limited to the leading 
order in the coupling strength A. However, potential 
applications of particle detectors will undoubtedly require 
calculations to higher orders in perturbation theory: some 
phenomena that have been investigated for quantized real 
fields using the original UDW detector, for example, rely 
on higher-order effects order (see Refs. [23, 46] among 
others). In preparation of investigations at higher orders, 
it will be convenient to derive a set of Feynman rules 
for each detector model. We will adhere to the following 
roadmap: In this section we develop the computation 
methods necessary to obtain Feynman rules. In Sec. VH, 
we will apply these methods to calculations up to second 
order in the transition amplitudes by way of example, and 
in Sec. VHI we will generalize the emerging pattern to 
Feynman rules. 

It is to be expected that new divergencies appear at 
each order in perturbation theory. As in any quantum 
field theory, there are three sources of divergencies in 
transition probabilities of UDW-type detector models: 

1. The amplitude of a single process can be infinite, as 
it was the case for the tadpolelike diagram in Fig. 1. 

2. Finite amplitudes can still amount to a divergent 
transition probability of the detector when the field 
is traced out: the sum over all the ingoing and 
outgoing field configurations can diverge. 

3. The perturbative series itself can diverge when tak¬ 
ing into account all orders. 

Notice, therefore, it is important to ascertain that those 
divergences can be regularized and renormalized away 
in order to be able to fully trust insights gained from 
using a particle detector model to describe the physics of 
a quantum field. 

Example: vacuum no-response probability. — In order 
to show which quantities we are going to need to evaluate, 
let us first consider a transition probability which involves 
the second order term in the Dyson expansion Eq. (25). 
Every application of the monopole operator p. Eq. (12) to 
one of the two energy eigenstates of the detector switches 
it back and forth between ground state \g) and the ex¬ 
cited state |e). In consequence, if the detector starts out 
in the ground state, only even orders in perturbation 
theory contribute to the probability for it to remain in 
the ground state, and only odd orders contribute to the 
probability for the detector to be excited. There are thus 
no corrections to the VEP coming from in Eqs. (25) 
and (26). Instead, we calculate the probability for the 
detector to remain in the ground state when interacting 
with a quantum field in the vacuum state, and call this 
probability the vacuum no-response probability (VNRP) 

Notice, however, that it would not be strictly necessary 
to go through the calculation in order to compute 
the VNRP: Since the same order perturbative corrections 
to the density matrix are traceless [5], it follows that 
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Po,g^g = 1 — Po,g^e- We will nevertheless compute it 
using to illustrate higher order techniques. 

Let F be the probed quantum field, and pitg) = 
10 ,( 7 ) ( 0 ,( 7 ! the density matrix of the coupled system at 
starting time tg. The VNRP is the ( 77 , g) component after 
tracing out the field: 


Po,g^g(tj to) — 

{g\tTF p{t)\g) = \{a,g]t\U{t,to)\0,g;to)\‘^ (106) 

a 


at time t > tg. In the limit to — 00 , t —)■ 00 , the actual 
duration of the interaction between detector and field 
is determined by the switching function y. Using the 
perturbative expansion Eqs. (25) and (26), the VNRP 
can be written as 


P 


o,g^git, to) — 1 + 2 ^ I^Re + 0{X^) , 

(107) 

where the transition amplitudes at order n are 

^(’;^ = (-iA)"(a,(7;t|C/(")(t,to)|0,ff;to) , (108) 


The contribution at order zero is only nonzero for the 
final state | 0 ,( 7 ) since 17° = 1 : 


Higher orders. — Generalizing the above example to 
order n and arbitrary initial and final states, we need to 
evaluate expressions of the form 

{g\[a-nTp{F)---pit^)][a+f\g) (114) 

for the detector, where a, & G {0,1} to allow for initial 
and final states jg) and |e). Similarly, for model 1 and 
quantized real fields we need to calculate 

(0|(ap)p[T<I>(27i)...cf>(y„)](4)fc|0) , (115) 

where (afc)fc = a\,_^ ... is an arbitrary product of cre¬ 
ation operators [such as, e.g., (al)k = 

(®I)fe ~ allowing for arbitrary number eigenstates 
content in the initial state, and (ap)p are annihilation 
operators for arbitrary final states. In case of model 3 
(for quantized complex fields), we need 

(0| (ap)p (bp)p [T : <I>l(?7i)$(j7i) : • • • 

• • •: ‘^Hyn)Hyn) ■■ ] (al)fe (&^)fc |0) (116) 

and in the case of model 2 (for quantized real fields), 
the previous expression simplifies because the field is self- 
adjoint, $1 = $, and there is only one type of ladder 
operators. Finally, for model 4 coupling to quantized 
spinor fields expressions of the form 


■^a.,g 


1 for a = 0 
0 else 


(109) 


As we have argued earlier, order n = 1 does not contribute 
at all, a2^1 = 0. Thus, only the amplitude at order n = 2, 


Ai^j = (-i)^ [ dti ( At2 

Jto Jto 

(a,5;f|Ti7i„t(ti)iJi„t(t2)|0,g;to) , (HO) 

remains to be evaluated to obtain the VNRP. 

As an example, consider the renormalized detector 
model 3 for quantized complex fields with interaction 
Hamiltonian Eq. (85). The second order correction is 


X f dti f dt 2 x{ti)x{t 2 ) {g\Tp{ti)p{t 2 )\g) 

Jto Jto 

X dyi dy2p{xo,yi)p{xo,y2) 

X (0|T: $'f(?7i)4>(j7i) :: $'f(y2)^’(272): |0) . (Ill) 

The main difficulty now rests in evaluating the time- 
ordered expectation values 

{g\Tp.{ti)p.{t 2 )\g) (112) 


(0| (ap,r)p.r {bp,f)p,f [T : 4'(77i)T(77i) : • • • 

• • •: T(27„)T(j 7„) : ] (a^ |0) (117) 

appear. 

In short, we need to evaluate the vacuum expecta¬ 
tion value of time-ordered products of field operators, 
which may contain normal-ordered subproducts, with a 
prepended product of annihilation operators and an ap¬ 
pended product of creation operators. Since the field 
operators and the monopole operator are essentially sums 
of ladder operators as well, a good strategy is to normal- 
order the entire product by successively commuting (or 
anticommuting) ladder operators. This procedure will 
furnish us with variants of Wiek’s theorem. Moreover, 
the presence of the time-ordering symbol leads to the 
appearance of Feynman propagators instead of ordinary 
commutators. This approach is closely related to standard 
techniques in elementary particle physics, where elements 
of the scattering matrix S are computed to obtain prob¬ 
abilistic predictions about the outcome of scattering ex¬ 
periments. There, the Lehmann-Symanzik-Zimmermann 
reduction formula similarly allows one to express elements 
of S in terms of vacuum expectation values of products 
of time-ordered fields [43, 47]. 

Time- and normal-ordering .— Let us briefly revise 
the definition of time-ordering and normal-ordering for 
bosonic and fermionic fields. In the bosonic case, the 
time-ordering symbol is defined as 


and 


(0|r:<l>t(j;i)(O(yi)::$t(y2)cI>(272):|O) . 


T$(a:)$(y) 

= - 27 °)$(a;)$(y) -k 6 ( 77 ° - a:°)$(y)$(a;) . (118) 
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In particular, exchanging two operators inside a time- 
ordered expression does not engender any change: 

T^{x)^{y) = r$(y)$(a;) . (119) 

Normal-ordering is achieved by moving all annihilation 
operators to the right: 

: Ofeo],: = ajjOfc , : a^ak : = a^ak ■ (120) 

For fermionic fields, however, 

= 0(x° - 2/°)vI-(x)vI/(y) - 0(yO _ x°)vl/(y)vl-(x) (121) 

so exchanging two fields introduces an additional sign 


9(cj) 



K(w) 


T'i/{x)'S{y) = -T'S{y)'i>{x) . (122) 

Similarly, a sign is introduced when exchanging ladder 
operators in normal-ordering: 

■ ^k,s^p^r ■ — • ( 123 ) 

A. The monopole moment as a field 

A two-level particle detector is, in a way, a fermionic 
system: the ladder operators Eq. (6) of the monopole 
operator satisfy the anticommutation relations 

{cr+,cr+} = 0, {cr+,cr"} = l, {cr",cr"} = 0. 

(124) 



(b) Contour in eq. (131). 


K(w) 


By comparison with the ladder operators of the quantized 
spinor field it is easy to see that, formally, a~, corre¬ 
spond to the ladder operators of a fermionic field which 
(1) has only one mode, (2) does not have spin degrees 
of freedom, and (3) has particle modes which are their 
own antiparticles. The analogy extends to the monopole 
operator: in the Heisenberg picture, 

y{t) = cr+e+‘^‘ -b , (125) 

which can be seen as a simplification of the mode ex¬ 
pansion of the quantized spinor field Eq. (A58). In the 
following, we will therefore regard the monopole operator 
as a fermionic quantum field ij.{x) = fx(t,x) = y{t) on 
(1, n)-dimensional Minkowski spacetime which is constant 
in space, rather than as an observable of a system in first 
quantization. The time-evolution of this held is described 
by a Klein-Gordon equation: from Eq. (125) immediately 
follows -b n^)/i(f) = 0, which can be rewritten as 

(□ + = 0 (126) 

since A^(a;) = 0. While this held is dehned on the entire 
spacetime, it interacts with the quantum held only for 
a limited time and in a certain place, according to the 
detector’s spacetime prohle /. 


Figure 3. Curves for the complex contour integrations in 
Eqs. (130) and (131). 


In terms of particle number eigenstates, this held is very 
simple: since there is only one mode, no spin quantum 
numbers and no distinct antiparticles, there can at most 
be one quantum; excitation of more quanta is forbidden 
by the Pauli exclusion principle. The two particle number 
eigenstates are simply the two energy eigenstates of the 
detector: the ground state |g) for no quantum, the excited 
state |e) for one quantum. The detector behaves as a sort 
of completely delocalized Grassman scalar [31, 48, 49]. 


B. Feynman propagators 

Deriving the Eeynman propagator Dp oi the UDW- 
type detector (or, more precisely, its monopole moment) 
runs analogously to the derivation of the Feynman propa¬ 
gator Sp of the quantized spinor held. We will treat the 
detector hrst, and generalize to the well-known Feynman 
propagators of the other helds for the sake of complete¬ 
ness. 
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1. UDW-type detector 


The Feynman propagator of the UDW-type detector is 
defined as the time-ordered vacuum expectation value 

DF{t,s) = {g\T^l{t)^l{s)\g) , (127) 

the resulting two 2-point correlators are 

{ p——s) f '> H 

for t < s ■ 


Note that Df{0) is ill defined, since Df —t ±1, depending 
on whether t < s or t > s. For t ^ s, we may simply 
express the Feynman propagator using the step function 
0 as 

DF(t,s) = . (129) 


This form of the Feynman propagator has the disadvan¬ 
tage that the exponential factor switches signs depending 
on which time is larger—in calculations of transition prob¬ 
abilities, one would therefore not end up with the Fourier 
transform of the spacetime profile / when integrating over 
all times in the Feynman rules later on. It is therefore 
useful to rewrite the expression using a complex contour 
integral (as it is usually for the Feynman propagator of a 
quantum field [47]). The two cases of time-ordering are 
then implemented as two different contours. 

By the residue theorem, we can rewrite 


g-in(t-s) 



D)(w -I- ft) 


(130) 


where the curve Fi is given in the first part of Fig. 3 and 
encloses the pole at w = +^1. Similarly, 


g-i-in(t-s) 



dio 

27r (w 


_ _ 

n){u} -|- D) 


(131) 


for r2 as shown in the second part of Fig. 3, enclosing the 
pole at u! = —O. In order to be able to integrate along 
the real axis, we shift the poles off the axis by ±ie and 
straighten out the curves as shown in Fig. 4. After the 
evaluation of the integrals, we will take the limit e —?► 0. 
Since r( = r(r) U ri(r) and F^ = r(r) U r2(r) as drawn 
in Fig. 4, the integral over each of the entire closed curves 
can be split into two integrals. In the limit r —)■ oo, the 
first integral simply turns into an integral over the entire 
real line. The second one, on the other hand, evaluates 
to zero in both cases: 



K(a;) 


Figure 4. Curves for the complex contour integration in 
Eqs. (132) and (133). 


By introducing time derivatives 

— in(t—s) _ iQ(t—s) +iQ(t—s) 


^ -i-in(t-s) 

n 


(134) 


and evaluating them inside the integral, we obtain 


DF{t — s) = 


lim 2 
£->■0 



duj 

2tt {lu 


lUJ 




D -f ie)(cu -f D — ie) 


(135) 


Since in the limit 


lim(a; — D -|- ie) (w -I- D — ie) = lim oj^ — +i6 , (136) 

£->■0 (S-J-O 

the above expression is abridged in the usual way to 


Dp{t — s) = lim 2 


duj 


lUJ 




£->0 J_ac, Stt oj^ — + ie 

(137) 

If the arguments of the propagator are exchanged, a 
sign is picked up. 


Df{s — t) = —DF{t — s) 


(138) 


0 (t — s) e = 0(t — s) lim 217 

£->■0 

r°° doj _ i_ , 

2 tt (w — D -I- ie)(a; -I- D — ie) 


and 

0(s - t) = 0(s - t) lim 2D 

£->■0 

r°° diu i 

7_oo 27r (w — D -I- ie) (w -I- D — ie) ^ 


(132) 


(133) 


since Tg{t)fi{s) = —Tfi{s)iJ,{t). 


2. Quantized scalar field 

Although well known in the literature, we repeat the 
definition of the Feynman propagator for quantized scalar 
fields: 


GF{x-y) = {0\T^x)^Hy)\0) . (139) 
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It can be expressed as 


1. Quantized complex field 


Gf{x -y) = 


lim 
€—>^0 Xj 


-E 


dk^ 


L„o 27r + : 




see, for example, Ref. [47]. Note that for x = y, 


Gf{o) = ^y — 
2 L" V Wfe 


00 . 


(141) 


Unlike for the detector, no sign is picked up under ex¬ 
change of the arguments: 


Gf{x -y) = GFiy - x) . 


(142) 


3. Quantized spinor field 


The Feynman propagator of the quantized spinor field 
is, by definition, 

- y)^ = (0|r4/-4(a;)TB(2/)|0) , (143) 

and can be expressed as [47] 


SF{x-y)B = 
1 


lim 

L’ 


E 


dk^ i[k • j + m 


2tt -I- : 


-ik{x-y) 


(144) 


The above expression is valid in (l,n) dimensions, for 
both n = 1 and n = 3; the massless case is correctly 
recovered by setting m = 0. 

It is straightforward to check that the Feynman propa¬ 
gator diverges in the coincidence limit 

S'i?(0) —>■ oo . (145) 


Unlike for the detector, where the exchange of the propa¬ 
gator’s argument only produced a sign, there is no simple 
relation between Sf{x — y) and SFiy — x): 


SF{y-x)^ 

= _ lini ^ V r ^ „-ik{x-y) 

e^o L" 4" J-oo 2^ k'^-m'^ + ie 

(146) 


C. Wick’s theorem 

In this section we will first address Wick’s theorem 
for quantized scalar fields. Then, we will formulate the 
fermionic version of the theorem for spinor fields, and 
subsequently simplify it in order to arrive to the case of 
the detector (monopole field). Notice that full proofs of 
the theorems are given in Appendix C. 


Ultimately, we want to calculate expectation values like 
Eq. (116), or more generally: 


(0| {ap)p {bp)p [T$ “^(a;!) • • • 

•••4>^"(a:n)](al)fc(4)fc|0) , (147) 

where the index indicates that the fields may be con¬ 
jugated or not. Since the vacuum expectation value of 
normal-ordered products of operators vanishes, it is useful 
to rewrite the original sequence of operators as a sum of 
normal-ordered expressions. To achieve normal-ordering, 
the fields need to be expanded in terms of ladder opera¬ 
tors, and the ladder operators commuted with each other 
until normal-ordering is reached. The well-known Wick 
theorem allows to do so systematically. We will proceed in 
two steps in evaluating the above expression. In the first 
step, we present Wick theorem for time-ordered products 
of field operators. In the second step, we then extend 
this theorem and allow for more ladder operators to be 
prepended and appended. 

It is useful to decompose the held operator into a part 
containing the creation operators and a second one con¬ 
taining the annihilation operators: 

4>(a:) = 4>o(a;) -f 4>|[(a;) , (148) 


where 

^a(^) = ^ak , ^b{x) = (149) 

k k 

and (fk are the mode functions of the quantised scalar 
held. 

Time-ordered produet of quantized complex fields .— 
The above separation facilitates normal-ordering series of 
helds. If there are only a few held operators, this can still 
rapidly be done by hand. For example, 

$(a;)$'l'(y) = : $(a;)$'l'(?/): -H[4>o(a;), $1(2/)] • (150) 

It is straightforward to check that the commutator is iden¬ 
tical to a nontrivial two-point function of the quantized 
complex held: 

[$a(a:), $1(2/)] = 1 (0|$(x)$l(2/)|0) . (151) 

Thus, time-ordering of the helds introduces the Feynman 
propagator, such that for example: 

T$(a;)$l( 2 /) = : $(a;)$l( 2 /): +1 Gf^x — y) . (152) 

These results are readily generalized if more helds are in¬ 
cluded. To keep the notation light, Feynman propagators 
are customarily written as contractions 

<^{x)<^\y) = tGF{,x-y), $(a;)$( 2 /) = 0 , 

$l(a;)$(2/) = 1 0 ^( 3 : — 2/) ) $l(x)$l( 2 /) = 0 . (153) 
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Theorem VI. 1 (Wick’s theorem for time-ordered quan¬ 
tized complex fields). A product of n time-ordered field 
operators of the quantized complex field which may con¬ 
tain subproducts that are normal-ordered, can be rewritten 
as follows: 

= : (0:2) •••$"" (Xn): 

+ ^ ; (x,) • • • (xj) • • • (xn) ■■ 

single 

contractions 

-f ^ (x,) • • • (xj) • • • 

double 

contractions 

■■■^^Hxk)---^^^ixi): + --- . (154) 


If it is even, it can be expressed as the sum over all full 
contractions: 

(0|(ap)p (6p)p ^(a;i) • • • $ "(a;„) ] (afc)fc (^^)fc|0) 

= X! “pi ■ ■ ■ f •••$"'= (a^fc) •• • 

all full 
contractions 

. (159) 

Kn 

The contractions between fields that do not belong to the 
same normal-ordered subproducts are 

r I I I, 

$'*'(a;i) • • • ^(xj) = 1 Gpixi - xj) = ^{x^) ■ ■ ■ $'(a;^) , 

(160) 

the contractions between ladder operators 


Two contracted fields belonging to different normal-ordered 
subproducts are to be replaced as follows: 


a/c' * * Up — 1 — bf^ • • • , 


#(xj) • • • 4>(xj) = 1 Gpixi - Xj) 

I I , 

^{xfi ■ ■ ■ 4>T(xj) = 1 Gpixi - Xj) 

I- 1 (155) 

$(xj) • • • 4>(xj) = 0 

$l(a;i) • • • 4)1 (xj) = 0 . 

Contraction of any pair of fields belonging to the same 
normal-ordered subproduct vanish. 

For three fields, for example, 


and, finally, the mixed contractions 

afc---4>l(x) = tT*k{x) , ^(x)---al = l(pk(x) 

bk---^(x) = Itfil(x) , 4>1 (x) • • • 6 j[, = l(pk(x) ■ 

(162) 

All other contractions vanish, in particular contractions 
between fields belonging to the same normal-ordered sub¬ 
products. 

A proof of this theorem is given in Appendix C 1. 


T 4 >(x):$l(y)$( 0 ): 

= : <l)(x)<i)l(j/) 4 >( 2 ;): -|- : <i)(x)<l)l(j/) 4 )( 2 ;) : 

= : <i)(x)<i)l(j/) 4 >( 2 :): -|-1 Gf{x — j/) 4 >( 2 ;) . (156) 

This theorem is a canonical result in quantum field theory 
[43, 47]. 

Theorem VI. 2 (Wick’s theorem for full expecta¬ 
tion values of the quantized complex field). Con¬ 
sider a time-ordered product of quantized complex fields 
T4)*^) (xi) • • • 4)'^" (x„) that may contain normal-ordered 
subproducts. The indices Ci € {o, f} indicate whether the 
field operator is conjugated (ti = f, such that 4)*^) = 4)1 j or 
not (ci = o, such that 4)"^) = 4)J. Then prepend products 
of annihilation operators 

{o.p)p = Opi • ■ • flpM ) {bp)p = bp^ ... bp^ (157) 
and append products of creation operators 

(4)fc = al • ■ • aL . (4)fe = 

The vacuum expectation value of the resulting product of 
operators vanishes if the total number of operators is odd. 


2. Quantized real field 

Note that the field operator can be separated as 

4>(x) = 4>a(x) -f 4-1 (x) . (163) 

Unlike for the quantized complex field, the commutator 

[4),(x),4)(y)] = ( 0 | 4 )(x) 4 )( 2 /)| 0 ) (164) 

does not vanish. This changes the contractions in the 
equivalent of theorem VI. 1 for the quantized real field: 
contractions between 4)(x) and 4)(y) do not vanish for 
the quantized real field. Similarly, the contraction be¬ 
tween Ofc and 4)(x) does not vanish in the equivalent of 
theorem VI.2, since 

[ak,Hx)] = lipl{x) 

+ (165) 

[4>(x),aj[,] = tifkix) 

are now nonzero. Keeping this in mind, we can immedi¬ 
ately state: 

Theorem VI.3 (Wick’s theorem for full expectation val¬ 
ues of the quantized real field). Consider a time-ordered 
product of quantized real fields T4)(xi) • • • 4)(x„) that may 
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contain normal-ordered subproducts. Then prepend prod¬ 
ucts of annihilation operators 

{ap'jp = Up^ ... apj^ (166) 

and append products of creation operators 

(al)fe = all ■ ■ • “L • (167) 

The vacuum expectation value of the resulting product of 
operators vanishes if the total number of operators is odd. 
If it is even, it can be expressed as the sum over all full 
contractions: 

(0|(ap)p[r$(i,i)...cO(a;„)](al)fe|0) 

= ^ Op, • • • ^{Xi) ■ ■ ■ ■ ■ ■ ^{Xk) ■ ■ ■ ^{xi) • • • al„ • 

all full 
contractions 

(168) 

The contractions between fields that do not belong to the 
same normal-ordered subproduct are 

$(xi) • • • $(a;j) = 1 Gpixi - Xj) , (169) 

the contractions between ladder operators 

Ofc • • • d], = 1 4,p , (170) 

and the mixed contractions 

ak-■-^{x) = l(plix) ^{x)---al = tipk{x) . (171) 

All other contractions vanish, in particular contractions 
between fields belonging to the same normal-ordered sub¬ 
products. 

3. Quantized spinor field 

For spinor fields, we need to work out how to evaluate 
vacuum expectation values of the form 


Note that this convention varies slightly from Eqs. (148) 
and (149) used for quantized scalar fields. 

Time-ordered products of quantized spinor fields. — Two 
spinor fields are easily normal-ordered by hand, for exam¬ 
ple: 

5'(a:)4'(j/) = : 4'(a;)4'(y): -^{^^(a;), 4'o(y)} • (175) 

The anticommutators can be evaluated using Eq. (A36) 
for massive fields in (1,3) dimensions [and the analogous 
equations for massless fields, and in ( 1 , 1 ) dimensions 
respectively]. It is straightforward to see that 

{T^(x),^a,B(y)} = l(0|vI/^(a:)^B(y)|0) . (176) 

The Feynman propagator appears again if the fields are 
time-ordered. For example, 

T'S{x)^iy) = : ^{x)^{y): +1 Spix - y) , (177) 

where we have used that : 4'(j/)4>(a:): = — : 4'(a;)4'(?/):. 

The generalization of the above is given by the fermionic 
version of Wick’s theorem for spinor fields: 

Theorem VI.4 (Wick’s theorem for time-ordered quan¬ 
tized spinor fields). Consider a time-ordered product of 
spinor fields (xi) • • • 4^*^" (a;„) which may contain 

normal-ordered subproducts. The indices e, G {0)~} indi¬ 
cate whether the field operator is Dirac conjugated (ci =~, 
such that = TJ or not (ci = o, such that 
It can be normal-ordered as follows: 

(Xi)'!/’'^ {X2) • • • T*'" {Xn) 

= :4'''i(a;i)4'‘'2(a;2)---4'""(x„): 

+ ^ : 4''^! (xi) • • • (xj) • • • 4''^" (x„): 

single 

contractions 

+ ^ : 4''^! (xi) • • • T'^‘(x,)---T'=^ (xj) ■ ■ ■ 

double 

contractions 

+ . (178) 


(0| (ap,r)p.r {bp,r)p,r ' ' ' 

•••'I'^"(^n)](al,Jfc..(6l,)s.,|0) (172) 

where once more the index indicates whether the field is 
conjugated (T), or not (T). Let us again split the field in 
two parts, according to the ladder operators, by defining 

^a(a;) = ^ ak,s 

k.s 

(173) 

'^bix) = _(x) , 

k,s 


Contractions of any pair of fields belonging to the same 
normal-ordered subproduct vanish. Two contracted fields 
belonging to different normal-ordered subproducts are to 
be replaced as follows: 


I- 1 

'^AiiXi)'^^^{Xj) = -ISpiXj -Xi)^l 

4'^*(a;i)4'.4, (xj) = +lSFixi- Xj)f^ 
4'(x,)4'(a;j) = 0 

I-1 

T(x,)4'(a;j) = 0 . 


(179) 


such that 


'I>{x) = 4'a(a;) -I- Th)^) 
^(a;) = ^^(a;) -f ^^(x) . 


(174) 


If contracted operators are not adjacent, the term has to 
be reordered until they are, with every exchange of two 
operators introducing an overall factor of (— 1 ) to the 
term. 
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The last requirement that contracted operators need to 
be adjacent means for example that 


: 4'(xi)^(a:2)^(a:3)^(a^4): 

I-1 

= (- 1 )^ : T(x4)^(xi) ^(x2)^ (a:3): (180) 

Since this is again a canonical result [43, 47], we abstain 
from giving the proof here. 

Theorem VI.5 (Wick’s theorem for full expectation 
values of the quantized spinor field). Consider a time- 
ordered product of spinor fields (a;i) • • • 4'*^" (x„) like 
in theorem VI.4- Then prepend products of annihilation 
operators 


If contracted operators are not adjacent, the term has 
to he reordered until they are, with every exchange of 
two operators introducing an overall factor of (— 1 ) to 
the term. All other contractions vanish, in particular 
contractions between fields belonging to the same normal- 
ordered subproduct. 

A proof is given in Appendix C 2. 

4 . UDW-type detector 

Last but not least, we would like to compute vacuum 
expectation values of the from given in Eq. (114), 

{g\[a-nT^,{C)■■■n{t^)][a+f\g) (187) 


(Op.r)p,r — ^(p,r)i ■ ■ ■ ^(p,r)M 

{bp,f)p,f = ^(p,f)i • ■ • ^{p,f)jQ ; 

and append products of creation operators 


(181) 


As before, we decompose the monopole operator into a 
creation and an annihilation part, 

/x_ (t) = e-'^* a- , PL+ (t) = a+ , (188) 


~ ^(k,s)i ■ ■ • ^(fc.s)jv ■ 

The vacuum expectation value of the resulting product of 
operators vanishes if the total number of operators is odd. 
If it is even, it can be expressed as the sum over all full 
contractions: 

(0| {ap^r)p.r (bp^r)p.r [ T4/ ^(xi) • • • 

•••4''^"(a;„)](a;[, (&^ j)fc,5 | 0 ) 

= ^ipx)i • • • (xk) ■ ■ ■ 

all full 
contract. 

. (183) 


such that 

ttit) = V-{t) + V+it) ■ (189) 

Time-ordered products of monopole operators .— For 
two monopole operators, normal-ordering is swiftly ac¬ 
complished and yields 

tiiti)ti{t2) = : fj,{ti)fj,{t2) : -I {fi-ih), iJ.+ {t2)} , (190) 

where 

{M-(ti),M+(i2)} = le-‘^(*-*^) . (191) 

Now consider what happens if the product of monopole 
operators is time-ordered. For example, 

Tfj,{ti)pL{t2) = : p.{ti)p.{t2) : +1 Dpiti - ^ 2 ) , (192) 


The contractions between fields that do not belong to the 
same normal-ordered subproduct are 


'^A,{Xi)'^^^iXj) = -1 Sf{Xj - Xi)^{ 

I- 1 

'^'^'{Xi)'^A,iXj) =+lSF{Xi- , 

the contractions between ladder operators 


(184) 


and the mixed contractions 

I- 1 

ak,s'^A{x) - 

^^(a^)«Ls = lV'fc,^,+ (a;) 

bk,s'i’^ix) = l'ipk,.,-ix) 

I- 1 

^4i(a:)6^ = , __^(x) . 


(185) 


(186) 


using the definition of the Feynman propagator Eq. (129). 
Since the monopole operator /x(x) (as a held) has basically 
the same structure as the quantized spinor held, we can 
directly conclude that the generalization to n operators 
is given by 

Theorem VI.6 (Wick’s theorem for time-ordered 
monopole operators). A product of n time-ordered 
monopole operators p, can be rewritten as a sum of normal- 
ordered products as follows: 


Tp{ti) ■ ■ ■ p{tn) 

= :piti) ■ ■ ■ p{tn) : 

+ Yi ■ • • A(^i) • --hitj) ■ ■■hitn ): 

single 

contractions 

+ Y ■ ■ ■ ■ A(^i)-"/^( 4 )---/l(ii) • • • fiitl) ■ ■ ■ 

double 
contractions 

• • • p{tn) : + • • • ■ 


(193) 
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The contraction of adjacent operators is defined as 

ii{U)ii{tj) = 1 Dpifi - tj) . (194) 

If contracted operators are not adjacent, the term has to 
be reordered until they are, with every exchange of two 
operators introducing an overall factor of (— 1 ) to the 
term. 

Similar to the spinor field, this last remark implies, for 
example, 

: h{h) T{t2)T{h)T {U) ■ = -'■ h{h)T{h) T{h)T {U) '■ 


5. Full expectation value 

When prepending zero or one annihilation operator 
a~ and appending zero or one creation operator and 
taking the vacuum expectation value, more contractions 
appear: 

Theorem VI. 7 (Wick’s theorem for the full expectation 
value of the detector). Consider a product of n time- 
ordered monopole operators pi with a € { 0 , 1 } annihilation 
operators and b G {0,1} creation operators. Its vacuum 
expectation value is 

{g\[a-nTfx{C)---pc{tr.)][a+f\g) 

= X! ■ ■ ■ lf{tj) ■ ■ ■ Till) • • • cr+ . 

all full 
contractions 

(196) 

The contraction between two adjacent monopole fields is 
defined as 

pfti}pi{tj) = 1 Dpifi - tj) , (197) 

the contraction between ladder operator and monopole 
field 

a~pi{t) = (198) 

and finally the contraction between ladder operators 

(F^+ = 1 . (199) 

If contracted operators are not adjacent, the term has to 
be reordered until they are, with every exchange of two 
operators introducing an overall factor of (—1) to the 
term. 

The proof is presented in Appendix C 3. 


D. Diagrammatic computation 

A usual way to visualize the contractions between field 
operators in the above theorems, as well as the next step 
towards Feynman diagrams, is the use of line diagrams. 
Operators in a given product are represented by points 
(vertices), and contractions are visualized as lines between 
the vertices. The vacuum expectation value of a given 
product of operators can be obtained in a structured 
manner by drawing such a diagram for every nontrivial 
full contraction (meaning that there are no wholly uncon¬ 
nected vertices), assigning the corresponding contraction 
as a factor to each line, and then summing up the terms 
obtained from each diagram. 

To distinguish them from the fully fledged Feynman 
diagrams that we will develop in Sec. VIII, we will in the 
following sometimes refer to these diagrams as auxiliary 
diagrams. 


1. UDW-type detector 

The following set of rules can be applied as a help in 
evaluating expectation values of the form Eq. (187) using 
theorem VI.7: 

Auxiliary diagrams for the detector. 

1. If a creation operator is present, draw a vertex on 
the left end of the diagram. If an annihilation op¬ 
erator is present, draw a vertex on the right end of 
the diagram. These vertices will he referred to as 
external vertices (^ingoing on the left and outgoing 
on the right). 

2. In the middle, draw a horizontal series of vertices, 
one for each ti, starting with tn on the left, and label 
them accordingly. These vertices will be referred to 
as internal vertices. 

3. To account for the nontrivial contractions Eqs. (197) 
to (199), connect all vertices by lines such that each 
external and internal vertex has one line attached 
to it. 

4-. Draw one diagram for every way to connect all 
vertices using the rules in 3, such that no vertex is 
left completely unconnected. If this is impossible, 
the expectation value is zero. 

5. Associate the following factors to each line in a 
diagram: 

• 1 to a line connecting two external vertices 

• to a line connecting an ingoing vertex 
with an internal vertex ti 

• to a line connecting an outgoing vertex 
with an internal vertex U 

• Dp(ti — tj) to a line connecting two internal 
vertices U, tj; since Dp{ti — tj) = —Dp{tj—ti), 
and the overall sign is determined in other 
ways, the order is irrelevant. 

6. To obtain the correct overall sign of a diagram, fall 
back to the interaction picture and reorder the op- 
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g —int2 


g+iOti 






Figure 5. The three diagrams belonging to the example for 
the monopole operator. 


erators such that contracted pairs are adjacent, as 
deseribed in theorem VI. 7. 

7. The vaeuum expeetation value is obtained by sum¬ 
ming up the values assoeiated with each diagram. 

Example .— Assume the detector is excited both at the 
beginning and at the end. Including second order correc¬ 
tions to the amplitude, one encounters the expectation 
value 


E = {Q\a- [Tp{h)pi{t 2 )]cr+\Q) . (200) 

There are three full contractions which can be written as 
follows: 

E= (0| a~ [Tli{ti)li{t2) ] |0) 

+ (0| [T/ r(ti)/t(t2)]g + |0) 

+ (0| [T/ r(ti)^ (t2)]g+ |0) 

+ a~ pL{ti)p{t2) 

^ ^ (201) 

according to theorem VI.7. Up to the signs, the same is 
obtained from the three corresponding diagrams shown 
in Fig. 5. 

2. Linear coupling to quantized real fields 

To simplify the use of theorem VI.3 to evaluate the 
vacuum expectation values of the form 



Figure 6. The two diagrams belonging to the example for 
linear coupling to quantized real fields. 


1. On the left end of the diagram, draw vertices for 
every ereation operator, and label them by their 
momenta. Similarly, draw vertiees for every anni¬ 
hilation operator on the right-hand side and label 
them by their momenta. As before, these vertices 
will be referred to as external vertices ( ingoing on 
the left and outgoing on the right). 

2. In the middle, draw a horizontal series of vertiees, 
one for each yi, starting with yn on the left, and 
label them aecordingly. Once again, these are the 
internal vertices. 

3. To account for the nontrivial contractions Eqs. (169) 
to (171), connect all vertices by lines in the following 
way: 

• both external and internal vertices have one 
line attaehed 

• direct connections between ingoing vertices are 
forbidden; direct connections between outgoing 
vertices are forbidden 

4-. Draw one diagram for every way to connect all 
vertices using the rules in 3. If this is not possible, 
the expectation value is zero. 

5. Assoeiate the following factors to each line in a 
diagram: 

• Sk.p to a line eonneeting two external vertices 
k and p 

• TkiUi) to a line eonneeting an ingoing vertex 
k with an internal vertex yi 

• TkiVi) ® connecting an outgoing vertex 
k with an internal vertex yi 

• Gpiyi — Vj) to a line eonneeting two internal 
vertices yi, yj; since Gpixi — xj) = Gpixj — 
Xi), the order is irrelevant. 

6. The vaeuum expeetation value is obtained by sum¬ 
ming up the values associated with each diagram. 

Example .— Add two annihilation operators in front of 
the two time-ordered fields. The calculations can still be 
performed by hand in a reasonable time, yielding: 


(0|(ap)p[TcI>(2/i)---‘i>(y„)](aI)fc|0) , (202) 

the following equivalent diagrammatic method can be 
used: 

Auxiliary diagrams for quantized real fields, lin¬ 
ear coupling. 


(0|afcap[T$(yi)$(y2)]|0) 

= Afc(l/2) Ap(yi) + Tkivi) Tpiy2) ■ (203) 

Indeed, applying the above diagrammatic rules results in 
the two diagrams depicted in Fig. 6 and gives the same 
result. 



27 


3. Quadratic coupling to quantized real fields 

If the detector model 2 coupling quadratically to a 
quantized real field is used instead, expectation values of 
the following type are encountered: 


(0|(ap)p[T:$(yi)$(?/i):--- 

. (204) 

Since there are now normal-ordered subproducts of fields 
with the same argument, the above rules have to be 
adapted. 

Auxiliary diagrams for quantized real fields, 
quadratic coupling. 

The rules for quantized real fields and linear coupling still 
apply, with one exception in point 3: internal vertices 
have two lines attached; these are not interchangeable. 

Example. — We evaluate an expectation value similar 
to the example for linear coupling. There are four ways 
to contract two helds 'I)(j/i), and then two ways to 

contract the remaining fields with the ladder operators. 
In consequence there are eight terms in total, yielding: 

(Olafettp [T : <I>(yi)4>(2/i):: ^{y2)^{y2) :]|0) 

= -iGF{yi - 2 / 2 ) [Tkiv^) T*pi.yi) + Tk(.vi) T*p{y 2 )\ ■ 

(205) 

The corresponding diagrams are drawn in Fig. 7. The 
same result can be obtained more economically by drop¬ 
ping the distinction between the two lines attached to each 
internal vertex, leading to merely two distinct diagrams in 
Fig. 8. To compensate, we assign an additional symmetry 
factor of four to the line between the two vertices, which 
counts the different ways the fields ‘l>(?/i) and 4* ( 2 / 2 ) can 
be contracted while now resulting in the same diagram. 
In accordance with the convention in particle physics, we 
will in the following use symmetry factors rather than 
distinguishing different legs of the same vertex. 


f. Quadratic coupling to quantized eomplex fields 

Theorem VI.2 is useful to evaluate the vacuum expec¬ 
tation values of the form 


(0| (ap)p (bp)p [T : $^(2/i)$(2/i) :••• 

...:$t(y„)ci,(y„):](at)^(6t)_|o) , (206) 

which arise when calculating the response of the quadrati¬ 
cally coupling detector model 3. Here, the existence of two 
independent types of ladder operators (and accordingly 
the fact that the field is no longer self-adjoint) has to be 
taken into account. This is most conveniently achieved by 
giving the connection lines between vertices a direction, 
that is, by turning them into arrows. The following rules 
can be used to apply the theorem: 





Figure 7. The eight diagrams belonging to the example for 
quadratic coupling to quantized real fields. In each of the two 
lines, all four ways to connect the two internal vertices are 
realized. In the upper line, one field ‘i?( 2 / 2 ) is always contracted 
with Ofc, and one field ‘l?(i/i) with Up. In the lower line, the 
inverse case is accounted for. 



Figure 8. Diagrams belonging to the example for quadratic cou¬ 
pling to quantized real fields, reduced to two distinct diagrams 
by dropping the distinction between upwards and downwards 
lines of internal vertices. To compensate, an additional factor 
four is assigned to their contraction. 


Auxiliary diagrams for quantized complex fields, 
quadratic coupling. 

1. On the left end of the diagram, draw vertices for 
every creation operator and and label them 
by their momenta with ki, k 2 , where the bar is used 
to mark antiparticle operators. Similarly, draw ver¬ 
tices for every annihilation operator Up^ and bp^ on 
the right-hand side and label them by their momenta 
Pi and P 2 . 

2. In the middle, draw a horizontal series of vertices 
for each yi, starting with yn on the left, and label 
them accordingly. 

3. To account for the nontrivial contractions Eqs. (160) 
to (162), connect all vertices by arrows in the fol¬ 
lowing way: 

• vertices representing a^ and b are origins of 
one arrow each 

• vertices representing a and W are endpoints of 
one arrow each 

• vertices representing t/i are origin of one arrow, 
and endpoint of one arrow each 

• direct connections between ingoing vertices are 
forbidden; direct connections between outgoing 
vertices are forbidden 

• connecting an internal vertex with itself is for- 





29 


them accordingly. Also, indicate the appropriate 
spinor index Ai. 

3. To account for the nontrivial contractions Eqs. (184) 
to (186), connect all vertices by arrows in the fol¬ 
lowing way: 

• external vertices representing and b are ori¬ 
gins of one arrow each 

• external vertices representing a and b^ are end¬ 
points of one arrow each 

• internal vertices representing pi are origin of 
one arrow, and endpoint of one arrow each 

• direct connections between ingoing vertices are 
forbidden; direct connections between outgoing 
vertices are forbidden 

• connecting an internal vertex with itself is for¬ 
bidden 

• arrow directions must not collide 

4-. Draw one diagram for every way to connect all 
vertices using the rules in 3. If this is not possible, 
the expectation value is zero. 

5. Associate the following factors to each arrow in a 
diagram: 

• Sk.p Ss^r to an arrow connecting two external 
vertices {k,s), {p,r) (or their antiparticle pen¬ 
dants ) 

• 'f’k s+i^) arrow connecting an ingoing 

particle vertex (fc, s) with an internal vertex 
{x,A) 

• to an arrow connecting an ingo¬ 
ing antiparticle vertex (fc, s) with an internal 
vertex {x, A) 

• ‘f’pr + a{^) to an arrow connecting an outgoing 
particle vertex {p, r) with an internal vertex 
{x, A) 

• tppj._{x) to an arrow connecting an outgoing 
antiparticle vertex {p, r) with an internal ver¬ 
tex {x, A) 

• Spix-y)^ to an arrow going from an internal 
vertex (y, B) to an internal vertex {x, A) 

6. To determine the overall sign of a diagram, fall back 
to the interaction picture calculations and order the 
operators according to the contractions, as described 
in theorem VI. 5. 

7. The vacuum expectation value is obtained by sum¬ 
ming up the values associated with each diagram. 

Example 1 .— Consider the following term, which allows 
only one nontrivial contraction according to theorem VI. 5: 

(0|r:^(yi)4-(2/i)::^(y2)«'(2/2)|0) 

= (0| T : AAiyi)' ^^iyi) :: ^ b(?/2)^^(?/2) |0) 

=-tr Spiyi - y 2 ) SFiy 2 - yi) , (212) 

where the trace is over spinor indices. The same factors 
are obtained diagrammatically from Fig. 11. 


Spiyi - y2)i 

y2,B i i yi,A 

SF{y2 - yi)A 


Figure 11. Diagram belonging to Example 1 for quantized 
spinor fields. 



Figure 12. The two diagrams belonging to Example 2 for 
quantized spinor fields. 


Example 2. — Consider the term 

E = {0\ak.sbp,r [r:T(yi)T(j/i) :: T(y2)^'(2/2) : ] |0) . 

(213) 

According to theorem VI.5, there are two nontrivial ways 
to contract. By reordering and writing out the contrac¬ 
tions it can be cast into the form: 

E = {-lfak,.s'I’A{yi)'^'^{yi)'^B{y2)bp,r'It^{y2) 

+ (-l)®afc,s^'ij(2/2)4'^(2/2)4'A(yi)&p,r^''^(j/l) 

= -A'4}k,s,+{yi)SF{yi - y2)'(’p.r-{y2) 
-A'ilJk,s,+{y2)SF{y2-yi)tpp,r-{yi) ■ (214) 

The two corresponding diagrams are shown in Fig. 12. 


VII. SECOND ORDER: VACUUM 
NO-RESPONSE PROBABILITY 

We now turn to calculating the VNRP Eq. (107) for all 
four detector models up to second order as an example. 
The emerging patterns can then easily be generalized 
to full Feynman rules. We have already determined the 
transition amplitudes in Eq. (109), and A^l vanishes. 
Only the second order amplitude Aa,g given in Eq. (110) 
remains to be evaluated. 



A. Linear coupling to quantized real fields 
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To calculate the VNRP for the original UDW detector, 
model 1 in Table I, we insert the interaction Hamilto¬ 
nian Eq. (14) in the amplitude Eq. (110). Combining 
the switching function and the spatial profile into the 
spacetime profile / gives the amplitude 


Dpiyi - 1 / 2 ) 
y 2 yi 



Gpiyi - 2 / 2 ) 


dyifiyi) J dy 2 !{y 2 ) 

X (5|r/i(yi)M(y2)|ff) ((('|r$(yi)«'(j/2)|0) , (215) 

where each integration runs over [to,t] x R", and \(f)) are 
energy eigenstates with energy uj^. Theorem VI.3 and 
theorem VI.7, and their diagrammatic representations 
can now be applied to evaluate these amplitudes for any 
given final state (p of the quantized real field. If the final 
state is the vacuum, Ip) = |0), there is only one way to 
contract all fields: 


= J dyifiyi) j dy 2 f{y 2 ) 

X ( 5 I Tli{yi)li{y 2 ) \g) (0| T4>(?/i)$(i/2) |0) . (216) 

The corresponding auxiliary diagrams for detector and 
field share the same vertices, and can be combined as 
shown in Fig. 13. Writing the propagators out introduces 
factors of the form —where d = (w, 0 ) —in case of 

the detector, and in case of the field. This allows 

to perform the integration over the yi] for example: 

J dj/i/(2/i)e-‘("+'=)*'i = f{d + k) , (217) 

where the domain of integration is [to, t] x R”. In the limit 
to —)■ — 00 , t —)■ 00 , this is just the Fourier transform of / 
as defined in Eq. (53). For simplicity, the same symbol 
/ will be used both for finite times and in the limit ± 00 . 
The amplitude then reads 


2iw 


doj 1 
2 ti L” 


E 


dk^ 

27r 


{-iX)f{d+k) 


a;2 — 112 _|- jg f^2 _ jjj2 _|_ 


(-iA)/(-d-fc) . (218) 


The factors and integrals that now additionally appear 
in the amplitude can be incorporated directly into the 
diagrammatic rules to obtain Feynman rules: For example, 
the last amplitude as a whole can be represented by 
the Feynman diagram shown in Fig. 14. To regain the 
amplitude from the diagram, we simply have to integrate 
(sum) the product of all the factors over the energies and 
momenta d, k running in the loop (with the appropriate 
measure). 

If the final state contains one particle, \(f>) = jlfc), the 
amplitude vanishes. 


- (0| akT^y,My2) |0) = 0 , (219) 


Figure 13. Linear coupling to quantized real field: auxiliary 
diagram for detector and field if final state \<j)) — |0). 


2iui 



{-iX)f{d + k) 


_ i _ 


Figure 14. Linear coupling to quantized real field: Feynman 
diagram at second order for the VNRP amplitude, final state 
\<P) = |0). 


since the presence of three operators does not allow for a 
full contraction. 

If the final state is \(j)) = jlfc. Ip), there are two ways 
to contract the operators in the field part: 


= (-iA)2e-H“'=+--)‘ 


X J dyifiyi) J dy2fiy2) {9\Tf4^y.{y2)\9) 


X 


(0| afcapT$(i/i)$(y2) |0)-h (0| afcapT$(yi)4>(?/2) |0) . 

( 220 ) 


We have already calculated these contractions in Eq. (203); 
the corresponding diagrams were drawn in Fig. 6. The 
resulting amplitude is 




1 


—iu>kt 




y/2uJkL^ 
“ duj 2iw 


A/2wpL" 

{-iX)f{d-p) {-iX)f{-d-k) 


l_^ 27r — ^7^ + ie 

+ (-iA)/(d - k) {-iX)f{-d - p)] , (221) 


where we have suppressed the limit e —)■ 0, the two corre¬ 
sponding Feynman diagrams are drawn in Fig. 15. Note 
that both diagrams essentially have the same amplitude 
since the two outgoing momenta k, p appear symmetri¬ 
cally. 

Final states \(f>) with more outgoing field quanta cannot 
contribute, because there are only two fields available to 
contract annihilation operators (outgoing states) with. 
Equivalently, on the level of Feynman diagrams, there 
are no further diagrams at second order. The vacuum 
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Figure 15. Linear coupling to quantized real field: Feynman 
diagrams at second order for the VNRP amplitude, final state 
\(j)) = |lfc, Ip). The vertex factors and the detector propagator 
are not shown. 


Figure 17. Quadratic coupling to quantized real field: two out 
of eight Feynman diagrams contributing to g. 


_ i _ 

—m'^ -\-iv 



{-i\)f{d + k+p) 


Figure 16. Quadratic coupling to quantized real field: Feyn¬ 
man diagram contributing to 


no-response probability is now obtained by inserting all 
amplitudes in Eq. (107); we refrain from writing it out 
here. 


B. Quadratic coupling to quantized scalar fields 


Recall that there are eight ways to fully contract this 
expression; the two corresponding diagrams are drawn in 
Fig. 8. The complete Feynman diagrams are drawn in 
Fig. 17 and have symmetry factor four. The sum of the 
corresponding amplitudes is 


A' 


( 2 ) 

t fc dp iP 

4 


g-iwfci g-icjpt 


duj 1 

y/2u}kL^ y'2ujpL^ J-oo 27r T" 
{-iX)f{d- p + q) 


E 

q ' 
2iu) 


— 2it 


— m? + id 
+ i-'A)f{d -k + q) 
i 


o;^ — + ie 

{-i\)f{-d-k - q) 
2iuj 

— Q2 -I- ie 


q2 _ 


-(-iX)f(-d-p-q) . (224) 


1. Quantized real fields 


For model 2 with the interaction Hamiltonian Eq. (85), 
there are three nonzero amplitudes of the form 

= (-iA)^e"“^-^ J dyifiyi) J dy 2 /(y 2 ) 

X {9\Tpiyi)p-{y2)\g) 

X {fi\T: HyiMyi) :: 4>(y2)$(2/2) : |0) . (222) 


Finally, if the final state has four particles, there are 
four ladder operators for the outgoing states that must be 
contracted with four field operators, which gives 4! = 24 
possible combinations. Figure 18 shows the six corre¬ 
sponding distinct Feynman diagrams. The symmetry 
factor for each diagram is four, and their amplitude totals 
to 


For final state |0), there are two ways to contract all field 
operators. There is a single corresponding Feynman dia¬ 
gram, drawn in Fig. 16, its amplitude therefore contains 
a symmetry factor of two: 


d — 2 

•^0,3 - ^ 


duj 1 V "v / 1 V "v I up 

Ott T.u 2-^ / 9-77- T.n 2-^ / P-TT 


dfcO 1 




00 27 r T 

k 


27r L” 


dp’ 

27r 


X (-iA)/(d + k+p) 


2iuj 


— Q2 _|_ ig p2 _ jjP _j_ 


p2 — r7j2 _|_ 112 


^i-iX)fi-d-k-p) . (223) 


If the final state contains two field quanta, the expec¬ 
tation value in the field part corresponds to Eq. (205). 




= 4 


g-iwfct g-iwpt g-iwqt g-iwjt 

V2WfcT” y/2ujpL^ y/2ujqL^ V2wjT" 


X 



dco 

27T 


{-iXf 




2iuj 

— -I- ie 


[Rd -p-q) f{-d -k-l) + Rd-p-l) R-d - k - q) 

+ Rd-q-l) f{-d -k-p) + f{d-k-q) R-d -p-l) 

+ f{d-k-p)f{-d-q-l) + f{d-k-l) R-d-p-q) . 

(225) 
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Figure 20. Quadratic coupling to quantized complex field: the 
two Feynman diagrams contributing to ^ 

J-fc) J-p ,9 


Figure 18. Quadratic coupling to quantized real field: six out 
24 Feynman diagrams contributing to iq ij s' 



Figure 19. Quadratic coupling to quantized complex field: 
Feynman diagram corresponding to 


2. Quantized complex fields 



Figure 21. Quadratic coupling to quantized complex field: the 
four Feynman diagrams contributing to j ^ 


When considering model 3 for the quantized complex 
held, there are again three nonzero amplitudes of the form 


y dyi/(yi) J dy2fiy2) 

X {9\Ty-iyi)fJ-{y2)\g) 

X . (226) 


If the hnal state is the vacuum, there is only one ■way to 
fully contract, unlike for the quantized real held. The 
expression in the held part was already evaluated in 
Eq. (207), the corresponding auxiliary diagram is given in 
Fig. 9. In combination with the diagram of the detector’s 
Feynman propagator, it becomes the Feynman diagram 
Fig. 19. Its amplitude is simply half of the amplitude in 
the case of a quantized real held as given in Eq. (223). 

If the hnal state of the quantized complex held contains 
one particle and one antiparticle, the amplitude admits 
two different ways to fully contract. The expression in 
the held has already been calculated in Eq. (210); the two 
corresponding Feynman diagrams based on Fig. 10 are 
given in Fig. 20. The resulting amplitude is 


A‘'^'> 

ifc.ip.g 


1 

4 


A 


( 2 ) 

Ife ? Ip 


(227) 


( 2^ 

differing from AlJ g in Eq. (224) only by lack of the 
symmetry factor. 


If the hnal state has two particles and two antiparticles, 
there are only four full contractions possible, as opposed 
to the quantized real held with 24 contractions. They 
correspond to the diagrams in Fig. 21, and the resulting 
amplitude is 


A' 


5 -icjfct 








( 2 ) ,_ 

V^uJkL^ y/2u}pL^ ^2a;qL" 

2ioj 


-(-wA- 
2 ti + ie L' 


fid-p-q) fi-d-k-l) 


+ f{d-q-l) f{-d - k - p) + f{d - k - p) f{-d -q-l) 
+ f{d-k-l)f{-d-p-q)\. (228) 


C. Quadratic coupling to quantized spinor fields 

For the spinor field detector, model 4, the missing 
amplitudes are of the form 

<g = (-iA)2e--^‘y dyi/(yi) J dy2fiy2) 

X {g\Tpiyi)p{y2)\g) 

X (V^IT: ^(j/i)^(2/i) :: ^(j/2)^(2/2) : |0) . (229) 

If the final state is |0), there is only one full contraction. 
The expression involving the field operators has already 
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k^ — m^+iS 


Figure 22. Quadratic coupling to quantized spinor field: Feyn¬ 
man diagram with factors corresponding to Note that 

the sign of the momentum in the spinor field propagator is 
such that it flows into the direction of the field line (fermion 
flow). 


been evaluated in Eq. (212) and yields twice the Feynman 
propagator of the spinor field. Writing it out according 
to Eq. (144), the amplitude can be expressed as 


- - 
”^0.3 - 


duj 1 


f~ dk° 1 


iiLU 1 V ^ / ui% r V ^ / dp^ 

O.ir JP' • ^ I O.ir TP • ^ I 977 


27r L 


27t 


27T 


X (—iA)/(c? — k + p) X 


2ia; 


i[f^ + m]^ 


i[p + m]^ 


+ w 


— Q? -\- ie k'^ — w? + 
{-'Pf{-d + k-p) , (230) 


where we have used Feynman slash notation f, = k ■ 
Combining the two auxiliary diagrams corresponding to 
this contraction yields the Feynman diagram Fig. 22. Note 
the overall sign —1 of the amplitude that arises from the 
reordering of the fermionic field operators. 

If there are two field quanta in the final state (particle 
and antiparticle), there are two ways to fully contract, 
analogous to the quantized complex field. The expression 
in the field part is similar to the one given in Eq. (214), but 
the annihilation operators are exchanged. This introduces 
an overall sign —1, such that 


A. 
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lfc,s 1 lp,ri5 
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ujkL 
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Tfl _ ^ A 


ujpL 
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«/ —OO 

2ia; 


duj 

27r 


^(-iA)/(d -k + q)— - „ . 

ZTT U)^ — le 


X —y^(-iA)/(-d-p-g) ■ (231) 

+ 10 



Figure 23. Quadratic coupling to quantized spinor field: the 
two Feynman diagrams contributing to 


The expression for massless fields is recovered by 
L~^ —)■ L“". The corresponding Feynman dia¬ 
grams are drawn in Fig. 23. 

The only other final state with nonzero probability 
amplitude contains two particle states and two antipar¬ 
ticle states. Note that the fermionic nature of the field 
requires (fc,s) {q,u) and (p, r) (Z,v) due to the 
Pauli exclusion principle. As for the quantized complex 
field in the comparable situation, there are four possible 
contractions—see Fig. 24 for the corresponding Feynman 
diagrams. For massive fields, the amplitude is 




(-iA)^ 
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0JqL^‘ 
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AiP- 




dio 


2 ia; 


27r P — + ie 


f{d-p- q) f{-d -k- l){Uq^n,AVpP{Uk,s,BVi^y) 
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+ f{d-q-l) f{-d -k- p){Uq^u,Avf^y){uk,s,BVp^r) 


+ f{d-k-p) f{-d -q- l)(uk,s,AVp^y)(uq^u,BVi^y) 
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p,r / 

A 


+ fid-k-l) f{-d -p- q){Uk,s,Apv)P,u,BVpP 

(232) 

and the case of a massless field is recovered like for 
Eq. (231). 


VIII. FEYNMAN RULES FOR DETECTOR 
MODELS 

Feynman diagrams and Feynman rules can be under¬ 
stood as an extension of the diagrammatic representa¬ 
tion of Wick’s theorem.In order to obtain a transition 
amplitude, there are additional factors and summations 
(integrations) to be performed. Since their occurrence 
also follows simple patterns, it is straightforward to in¬ 
clude them into the diagrammatic rules: Combining the 
auxiliary diagrams of the different field types into one 
diagram containing all field types gives us a Feynman 
diagram, and incorporating all additional operations into 
the rules for evaluating auxiliary diagrams gives Feynman 
rules. 
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Figure 24. Quadratic coupling to quantized spinor field: the 
four Feynman diagrams contributing to tit- 

Note that the formalism presented here does not make 
use of asymptotic free states, in contrast to the usual 
Feynman rules of particle physics. Rather, ingoing and 
outgoing states are tensor products of Fock states of the 
noninteracting systems. 

A. Towards Feynman rules 

Feynman diagrams .— Consider, for example, the ex¬ 
pression 

{g\ a~Tg.{yi)fj,{y 2 )a~^ \g) 

X (0| bp^rak.sT : :: ^(j/2)«'(2/2): |0) . (233) 


fields, and (massive) spinor fields from left to right: 

i 

a ;2 — 112 _i- ig ’ _ ^2 -(- ig ’ k'^ — m? + ie 

Each vertex contributes f{k) where k is the energy- 
momentum flowing into the vertex. Finally, all inter¬ 
nal momenta (the ones in the propagators) need to be 
integrated or summed according to 



for the detector (left) and the fields (right). 

Let us now consider the factors associated with incom¬ 
ing quanta on one hand, and cases where two ladder opera¬ 
tors are contracted, i.e., diagrammatically, when an incom¬ 
ing line is directly connected to an outgoing line. The con¬ 
tractions related to the first case are a~y.{ti) = 1 

I-1 

for the detector, Ofc • • • 4>l(x) = 1 Fki^) ^ quantized 
real field, and so on. The exponential functions simply 
contribute to the Fourier transform of /, such that the 
incoming momenta appear in the momentum balance 
of the vertex they are attached to. The normalization 
factors and spinors are left over and assigned to the in¬ 
coming diagram leg. Moreover, since we are working in 
the interaction picture, there is a time-dependent phase 
in addition to the ladder operator creating the state, for 
example, 

|e;to) |e), lU; <o) = |lfc) . (237) 

Altogether, the incoming legs of the diagram contribute 
factors 


The detector part allows for three different contractions, 
giving rise to three terms and the three corresponding 
diagrams drawn in Fig. 5. The field part has two possible 
contractions, with diagrams shown in Fig. 12. Since these 
two parts are multiplied, one arrives at six summands 
in total. Both the diagrams for the detector and the 
field always have the same amount of vertices (equal to 
the order in perturbation theory), and can be combined 
into six Feynman diagrams which simultaneously show 
contractions for the detector and the field. 

Feynman rules .— From the VNRP calculations in 
Sec. VII, we can directly deduce that outgoing quanta 
contribute the factors 
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(234) 
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(238) 


in the same order as before. 

If two ladder operators are contracted, they merely 
contribute Kronecker deltas ensuring that both are of the 
same type (momentum, spin) according to Eqs. (161), 
(170), (185) and (199). In the full amplitude, they addi¬ 
tionally contribute phases from the interaction picture. 
For example, in case of a quantized real field: 


{lk,g;t\U‘^°'>{t,to)\lp,9;to) = 4,p - (239) 

Similar results hold for the monopole moment, and the 
remaining fields. Disconnected lines from the left to the 
right end of the diagram will thus contribute 


to the amplitude [for the detector, quantized scalar fields, 
quantized spinor field particles and (massive) spinor field 
antiparticles from left to right and top to bottom]. In¬ 
ternal lines between vertices give (Fourier transformed) 
Feynman propagators as usual; for the detector, scalar 


^ ^ ^k,p 5 ^ ^k.p ^s,r 

(240) 

for the detector, quantized scalar fields, and quantized 
spinor fields, in that order. 
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The Feynman diagrams presented in this chapter, as 
well as the auxiliary diagrams in Sec. VID, do not contain 
lines for the detector in its ground state—in this they 
differ from the representations we used in Secs. Ill and IV 
and it differs from the diagrams used in [21]. The reason 
for this is that Feynman diagrams generally describe 
the energy-momentum flows involved in the interaction 
processes. The detectors (as macroscopic physical objects) 
do carry energy and momentum, but in the UDW-type 
detector models, only the excitations of its internal degree 
of freedom couple to the field. This is why lines are drawn 
not for the detector itself but instead for the excitation, 

i.e., for the quantum of energy, that the detector can 
carry. 


B. Feynman rules for single detectors 


We summarize the Feynman rules for all four UDW- 
type detector models considered in this work. The rules 
are applicable for detectors on Minkowski spacetime , the 
trajectory of the detector is determined by its spacetime 
profile /: if the spatial dependence varies with time in 
the restframe of the field, the detector is moving. Assume 
the field is initially in the number eigenstate \ai) at time 
to, and the detector in the state |/?i) (which may either 
be the ground state \g) or the excited state |e)). The aim 
is to calculate the probability amplitude 


{af,l3f,t\U{t,to)\a„l3i;to) (241) 


and add the resulting amplitudes to obtain A: 

n ^ 

A = J2 ( 243 ) 

k=0 

The four sections below specify how to draw and eval¬ 
uate the Feynman diagrams at a given order k for the 
different detector models and fields. 


1. Linear coupling to quantized real fields 


The Feynman diagrams for the detector model 1 cou¬ 
pling a single UDW-type detector to a quantized real 
field consist of the seven different types of elements listed 
below. Each element is associated with a factor in the am¬ 
plitude A[^'^ of an individual diagram. The corresponding 
factor is printed to the right of each element. Straight 
lines represent excitations of the quantized scalar field, 
wiggly lines excitations of the UDW-type detector. 


Elements in Feynman diagrams. 


2. k - p 4,p 


for the detector to end up in the state |/?/) and the 
field in |a/) at time t > to to order A" in perturbation 4- 

theory. In general, this is accomplished by combining the 
perturbative expansion Eqs. (25) and (26) of the time- 
evolution operator with the appropriate Wick theorem. 5- 
The calculations can be organized using a set of Feyn¬ 
man rules; the procedure itself is independent of the 
detector model and the type of field: 6. 



1 p-l-iwfcto 
y/2u)kL^ 


e 


— iQt 


1 

y/2uJkL'^ 


Generic Feynman rules. In order to ealculate a prob¬ 
ability amplitude of the type Eq. (241) to order A” in 
perturbation theory, proceed as follows: 

1. At a given order k < n, draw all possible Feynman 
diagrams using the rules appropriate for the detec¬ 
tor model and field in question, as described in the 
following sections. 

2. Evaluate the amplitude Al of each diagram by first 
multiplying the factors associated with the individual 
parts of that diagram, and then summing (integrat¬ 
ing) over all undetermined energies and momenta. 

3. Sum up the amplitudes at order k: their contribution 
to the total amplitude A is 

(242) 

i 

4-. Repeat the above steps for every order 0 < k < n. 


i-iX)fid-k) 

Here, 

f{d -k)= I dx° j da; f{x) (244) 

J tQ J IR,^ 

is the temporally cutoff Fourier transform of the detector’s 
spacetime profile /. Momentum flowing into a vertex has 
sign +I and is indicated by an arrow pointing towards 
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the vertex, momentum flowing out of a vertex has sign 
— 1 and is indicated by an arrow pointing away from the 
vertex. Depending on the situation, the direction of any 
momentum in the vertex (element 9) could of course be 
reversed. 

Constructing a Feynman diagram. To construct a 
Feynman diagram at order k and evaluate its amplitude 

(k) 

A\ , adhere to the following rules: 

1. Draw a horizontal series of k vertices. These will 
he referred to as internal vertices. 

2. If the detector is initially excited, \j5i) = \e), draw 
a wiggly line entering the diagram from the left 
(element 3). 

3. For every quantum in the initial state of the field 
\ai), draw a plain line entering the diagram from 
the left (element 4)- 

4- If the detector is excited in the end, \j5f) = \e), 
draw a wiggly line leaving the diagram on the right 
(element 5). 

5. For every quantum in the final state of the field |a/), 
draw a plain line leaving the diagram on the right 
(element 6). 

6. Connect the ingoing lines with outgoing lines (el¬ 
ements 1 and 2) or vertices, and the remaining 
outgoing lines with vertices as well. Connect the 
vertices with each other (elements 7 and 8). In the 
end, eaeh vertex has to have one deteetor line and 
one field line attached (element 9). 

7. If this is impossible, .4-^^ = 0 and = 0. 

8. Else, multiply the faetors eorresponding to the ele¬ 
ments in the diagram, integrate over the momenta 
running in the propagators (elements 7 and 8) by 


Constructing a Feynman diagram. Construetion of 
a Feynman diagram at a given order k proceeds as before, 
with the exception of 

6’. [.. . ] In the end, each vertex has to have one de¬ 
tector line and two field lines attached (element 
9’). 

and the appearance of symmetry factors because each ver¬ 
tex corresponds to two identical field operators ^{yi)^{yi): 

9. Multiply the amplitude with the symmetry factor of 
the diagram (compare Sec. VID 3). 


3. Quadratic coupling to quantized complex fields 


In case of the (leading order renormalized) model 3 
and a quantized complex field, one additionally has to 
distinguish between particle and antiparticle quanta. Dia- 
grammatically, this is achieved by turning held lines into 
arrows. 

Elements in Feynman diagrams. 


1 . 


g-in(t-io) 


2 . 





3. 



g+iHio 


I 




1 



(245) 



e 


— iQt 


and finally take the limit e —>■ 0 for any propagator g 
present to obtain the amplitude of the diagram. 




1 


2. Quadratic coupling to quantized real fields 


7. 



2iij 

w2_n2+ie 


The Feynman rules for the detector model 2 coupling 
quadratically to quantized real fields through Eq. (85) 
(renormalized at leading order) are largely equal to the 
ones given in the previous section. 


Elements in Feynman diagrams. 

Having quadratic coupling instead of linear cou¬ 
pling merely changes the vertex, all other elements 
are identical to the ones listed in the previous section: 



(-iA) f{d-k- p) 



Constructing a Feynman diagram. For quantized 
complex fields, step 6. from Sec. VIIIB1 has to be 
amended once more: 

6’. [.. . ] In the end, eaeh vertex has to have one de¬ 
tector line and two field lines attached (element 9). 
Arrow directions must not collide. 
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Unlike for the quantized real field, there is no need for in¬ 
troducing symmetry factors, since the two field operators 
underlying a vertex, are not identical (see 

Sec. VID 4). 


Elements in Feynman diagrams — massless field. 



1 „+iiA>htn 


U 


A 

k,s 


4- Quadratic coupling to quantized spinor fields 



1 

vW 




^k,s,A 


For a quantized spinor field coupled to an UDW-type 
detector through Eq. (84), the same elements appear, but 
are assigned different factors: 


Elements in Feynman diagrams — massive field. 9’. 

11 ’. 

2. -►— iP^r) 

k 



1 ^+iu:kto ^ 




1 




V 


A 

k,s 


k^+ie 


3. 



I 



5. 



6 . 





g-t-if2io 


ra g-t-iojfcfo 


m p+iujkto 

fiuJiQLA 


'^k.s,A 


Constructing a Feynman diagram. For spinor fields, 
the steps from Sec. VIIIB 3 are still applicable, but it is 
convenient keep track of spinor indices in addition: 

1 ’. Draw a horizontal series of k vertices, and label 
each one with a spinor index Ai. These will be 
referred to as internal vertices. 

Note that these rules do not provide the overall sign of 
the amplitude of an individual diagram at order k. 
As for the auxiliary diagrams in Sec. VID, the sign is 
most easily determined by writing out the underlying 
contractions and reordering the operators explicitly. 


7. 



e 


— iQ,t 


C. Feynman rules for multiple detectors 





’^k.s,A 


9. 



ni p+iujhtn ...A 


10 . 



2icj 

—n^-t-ie 


11 . 



i[^+ra\^ 

k^-m'^+ie 


When dealing with massless fields, the redefinition of 
the normalization factors in the mode functions (see Ap¬ 
pendix A 2) requires redefinition of part of the Feynman 
rules just stated: 


For questions in relativistic quantum information the¬ 
ory, often two or more detectors coupled to the same 
quantum field F are considered. Each detector i has its 
own Hilbert space, and a monopole operator acting 
on it. For N detectors, the Hamiltonian of the entire 
system is 


= + (246) 

(i) 

where each is the Hamiltonian of a free detector 
as given in Eq. (5), and each detector by itself interacts 
with the field as before. The above Eeynman rules readily 
generalize to such situations: each additional detector is 
equivalent to an additional monopole field pS’'\x). Be¬ 
cause each of these has its own Hilbert space, quanta (and 
the corresponding lines in Feynman diagrams) of different 
detectors always need to be distinguished. Since there 
are only two fields (for linear coupling) or three fields 
(for quadratic coupling) present in each interaction Hamil¬ 
tonian, there are still only two- or three-legged vertices, 
respectively. While direct interactions of the detectors 
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Figure 25. Example of a leading-order process involving two 
different detectors, blue and red (upper and lower wavy lines), 
quadratically interacting with a quantized complex or spinor 
field. Initially, the blue (upper) detector is excited. By ex¬ 
change of a field quantum, the excitation is swapped to the 
red (lower) detector under emission of a particle-antiparticle 
pair. 


do not occur in this model, they can of course interact 
indirectly through the quantum field they probe—Fig. 25 
demonstrates such a process. 


IX. CONCLUSIONS AND OUTLOOK 

A key feature of Unruh-DeWitt detectors, i.e., of first 
quantized systems that are coupled to quantum fields, 
is that they can be localized in space and time. UDW 
detectors therefore provide an excellent tool for extracting 
spatiotemporal information from quantum fields, at least 
in gedanken experiments. So far, most studies have used 
the simplest kind of Unruh-DeWitt detectors, namely 
detectors that couple linearly to a bosonic real scalar 
field. The use of those Unruh-DeWitt detector models has 
yielded important results, for example about the Unruh 
and Hawking effects [2, 15, 50-52], and, more recently, 
about the extraction of entanglement from the vacuum 
and communication via quantum fields [4, 5]. 

These successes motivate the development of Unruh- 
DeWitt detector models also for the detection of quanta 
other than real scalar bosons, such as fermionic and com¬ 
plex quanta. There are pioneering studies in this direction 
[17, 20, 32, 33], but so far they have not gone further than 
the calculation of transition rates. Indeed, we found that 
a direct application of these models beyond the computa¬ 
tion of transition rates leads to divergencies that cannot 
be regularized in the usual fashion, i.e., through spatial 
smearing and smooth switching functions. 

Here, we introduced improved Unruh-DeWitt type de¬ 
tector models that cure these divergencies by means of 
renormalization. The new UDW detector models therefore 
possess the full versatility of the usual bosonic Unruh- 
DeWitt detectors. For example, they now allow for the in¬ 
vestigation of the entanglement structure of the fermionic 
vacuum. 

To this end, we started by analyzing the four different 
particle detector models summarized in Table I. We then 


showed that, without further treatment, the vacuum ex¬ 
citation probability is already divergent at leading order 
in perturbation theory for the models 2-4, which are cou¬ 
pling quadratically in the field. These divergences cannot 
be regularized by the introduction of a spacetime profile 
(switching function and spatial smearing), unlike diver¬ 
gencies appearing in the case of the traditional linearly 
coupled UDW detector ]13, 20, 34, 35]. 

We found the origin of these divergences to be the fact 
that the quadratic coupling gives rise to loop diagrams (see 
left side of Fig. 1) which are the exact equivalent of tadpole 
diagrams in quantum electrodynamics (QED) ]45], with 
the detector playing the role of the electromagnetic field. 
Motivated by the analogy to QED we renormalized models 
2-4 at leading order in perturbation theory by normal¬ 
ordering the interaction Hamiltonian. Physically, this 
corresponds to setting the (otherwise infinite) interaction 
energy to zero. 

Finally, we derived Feynman rules for all four detector 
models in Sec. VHI. The Feynman rules allow for the 
systematic, graph-based calculation and renormalization 
of, for example, transition probabilities to arbitrary order, 
in the way familiar from particle physics. 

With the Feynman rules for the various detector mod¬ 
els at hand, one can now straightforwardly study the 
behaviors not only of bosonic but also of fermionic fields, 
for example, concerning the communication between two 
first quantized quantum systems (UDW detectors) via a 
quantum field, or the harvesting or farming of entangle¬ 
ment from the vacuum. In particular, the latter may now 
also allow one to develop quantitative tools to measure 
entanglement in the fermionic vacuum. 

There are various further promising possibilities, in 
particular for applying the new renormalized detector 
model for spinor fields. Among others, this detector 
model can now be used to study the fermionic Unruh 
effect beyond the detector transition rates. This should 
allow one to study the universality and the thermalization 
of a detector in the fermionic Unruh effect [46] and to 
study possible finite time effects (see for instance ]53]). 

We also addressed the question of which bosonic detec¬ 
tor model allows for a fair comparison with a fermionic 
model. We found that the leading order behavior of the 
renormalized models 3 and 4 can be compared on equal 
footing. 

Finally, while UDW detectors are of great value for 
gedanken experiments, let us recall that the UDW detec¬ 
tor for bosonic fields can be viewed as an idealization of a 
realistic system, such as the lowest two energy levels of an 
atom or molecule ]27, 28]. The atom detects photons by 
internal excitation that can later be measured via stan¬ 
dard experimental procedures. It should be interesting to 
explore to what extent the UDW detector for fermionic 
fields, model 4 in Table I, can also be viewed as a sim¬ 
ple idealization of a realistic system, such as the optical 
cavity envisioned in the introduction, which would detect 
fermions through electromagnetic mode excitations. More 
generally, it should be interesting to further investigate 
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how the proposed detector models along with their mode 
selection and switching functions could find experimental 
realizations. 


Appendix A: Notes on quantum fields 

In this section, we summarize results of the classical 
and quantum field theory of scalar fields on Minkowski 
spacetime, as well as spinor fields on and All 
fields are spatially restricted to the cavity 

[-L/2,i/2]" C R” (Al) 

with periodic boundary conditions 

^{t,x + S) = VdGL-Z" (A2) 

The objective is to establish notation and provide a refer¬ 
ence basis for the main body of this paper. These results 
are obtained in close analogy to the standard treatment 
of unrestricted fields, see e.g. [43, 47]. 


with respect to the Klein-Gordon product 

($ 1 , $ 2 ) = -i / x) [5o$2(t, x)] 

Jb 

— [9od>i(t, £c)] d>2(t, a:)da7 . (A9) 

It is useful to define purely spatial mode functions 

" . (AlO) 

These are orthonormal with respect to the T^-scalar prod¬ 
uct 

7 T^p) — / (^)v^p— ^k.p • 

Jb 

To avoid cluttered notation, we will sometimes drop 
the tilde distinguishing time-dependent from time- 
independent mode functions. 


b. Elements of canonical quantization 


1. Scalar fields 

a. Results in classical theory 

Consider a real or complex field $ on Minkowski space- 
time. The dynamics of this field are determined by the 
Lagrangian densities 

c% = 

(for the real and the complex field, respectively), which 
both yield Klein-Gordon equation 

=0 (A4) 

as equation of motion. A solution space to the Klein- 
Gordon equation is spanned by the family of mode func¬ 
tions 


‘Rk{t,x) 




(A5) 


where the momentum takes the following values: 


k e 271IL ■ Z” . (A6) 


Canonical quantization of the free real field leads to 
the mode expansion 

^(^7 x)='^ak ifkit, x) + al x) (A12) 

k 

for the Heisenberg picture field operator in terms of the 
annihilators and creators al- In the Schrodinger pic¬ 
ture, 

$(£c) = ^ 

\/2xJk 

The field operator is self-adjoined, 4)^ = $ , and the ladder 
operators satisfy the canonical commutation relations 

[afc,a^] = 14,p 
[ak,ap] = 0 = [aT,aT] . 

After renormalizing the vacuum energy by normal¬ 
ordering, the dynamics of the system are determined 
by the Hamilton operator 

= '^ujkalak ■ (A15) 

k 


ak Vk{x) -l-a][, ip*k{x] 


(A13) 


The energy ojk is determined by the dispersion relation 
ujk = i/lfcp -p m? . (A7) 

These mode functions have been normalized according to 
the condition 


(T^/c 7 ^p) — ^k,p 


(AS) 


Particle number eigenstates are created by acting on the 
vacuum state with creation operators 


l(nk)k) 


n 


1 


(«!)”'= io). 


(AI 6 ) 


Creation operators commute, so their ordering is not 
relevant. 
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Quantization of the free complex field leads to the mode 
expansion 




L 


dj^ c 


—iiOkt 


Vkix) 


+ (A17) 


in the Heisenberg picture. The field operator is not self- 
adjoint, the ladder operators satisfy the canonical com¬ 
mutation relations 


[ofe^Op] = 14,p , [bk,bl] = ISk^p , (A18) 


representation of Cliff(1,1) on C"* very similar to the 
Dirac representation. 

The dynamics of the field are governed by the Lagrange 
density 


= 'S{ij'^dp - m)'S , 

(A23) 

leading to the Dirac equation: 


— m)'S{t,x) =0 . 

(A24) 


Here 7 ^ are the Dirac matrices. Solutions to the Dirac 
equation are linear combinations of the mode functions 


and all other commutators not yet fixed by these relations 
vanish. The Hamilton operator of the free complex field 
is 

i? = ^ Wfc (alok + b^bk^ , (A19) 

k 


■^k,s,e(t,x) ='4}k,s,e(.x)e (A25) 

where s G {—V^} describes the spin degree of freedom, 
and k and e can take the same values as in Eq. (A 6 ). These 
mode functions are taken to be normalized according to 
the condition 


and particle number eigenstates are again generated by 
applying the creation operators on the vacuum. The 
complex field has 17(1) charge, so we need to distinguish 
between particle states created by a and antiparticle states 
created by b: 


\{‘^k)k: {j^p)p) — 


n 




( 4 rn 



(A20) 


2. Spinor fields 

a. Results in classical theory 

Let T be a spinor field on Minkowski spacetime 
or R^’^, taking values in C"* as a representation space 
of the Clifford algebras Cliff(l,l) or Cliff(l,3), respec¬ 
tively. In (1,3) dimensions we choose the irreducible Dirac 
representation, generated by the Dirac matrices 


(V^fc,s,e 5 V^p,r,(5) — '^fc.p ^€,(5 — (V^fc,s,e , ^p,r,(5) 5 

(A26) 

where we use the (generalized) scalar product 

{'Si, ^ 2 )= [ 'S\{t,x)'S2it,x)dx . (A27) 

JB 

Again, to lighten notation, we will sometimes drop the 
tilde if safely possible. 

The solution of the Dirac equation in the cavity B can 
be accomplished in close analogy to the case of a free 
spinor field, a detailed treatment of the latter can e.g. be 
found in [43]. The actual form and the properties of the 
mode functions depend on the spatial dimension as well 
as on the mass parameter m. 

(1,3) dimensions .— In case of a finite mass to 7 ^ 0, 
the mode functions in three spatial dimensions can be 
decomposed as 


tpk,s,€{x^ 



+eikx 

e , 


(A28) 


[12 

0 

1 

0 CTj 

' 0 
_1 

— 12 

II 

-Ui 0 


(A21) 


where II 2 is the two-dimensional identity matrix, and 


CTl = 





(A22) 


are the Pauli matrices. The Dirac matrices satisfy the 
anticommutation relations 


{7^7"} = 2r?^"l2 

for yL,v & {0,1,2, 3}, and therefore generate a represen¬ 
tation of the Clifford algebra Cliff(l, 3)® on C"*. In (1,1) 
dimensions, we use 7 ° and 7 ^ to generate a (reducible) 


where the spatial dependence has been separated from 
the spinorial part. For the latter, we will use the well- 
established notation Ufc,£,,e=+i = Uk,s and Uk,s,e=-i = 
Vk,s- Explicitly, they are 


tl'k.s — 


t'k.s — 


ujk + m 
2m 


LOk +m 
2m 






where cr* are the Pauli matrices, and 

= ( 0 ) ’ = ( 1 ) 


(A29) 


(A30) 
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are eigenvectors of a^. Multiplying two spinors yields 
1 


^k,s ^p,r — 


2m 


^ (wfc + m){ujp + m) 


2si{k^p^ — krp^) + k ■ p 


+ 


yj{ujk + rn){ijjp + m) 

1 2s{k}p^ — k^p^) — i{k‘^p^ — k^p^) 
2 w y/{ujk + m)(wp + m) 


5s,-r (A31) 


^k,s '^p,r — 


{(jJk + m){ujp + m) 

2m 


/ 2 sp^ 

2sk^ \ 

\Wp - 1 - m 

LOk + m) 

^ — 2sip‘^ 

k^ — 2sik‘^ 

ojp + m 

LOk + m 


5s.-r 


and 


'^k,s ^p,r — ^k,s ^p,r 


^k,s ’^p,r — {^k,s ’^p,r) 


(A32) 

(A33) 

(A34) 


The bar indicates the Dirac conjugate u = as usual. 
In particular, the normalization is chosen such that for 
fe = p all products simplify to 


^k,s,€ ^k,r,S — ^ 5s r 5s S • 


(A35) 


Moreover, products with the inverse order of spinors (such 
that they multiply to a 4 x 4 matrix) can be brought into 
a compact form as well: 


^k,sUk,s,B = ^ [^ ■ ^ + 

S 

^ts'^k,s,B = ^[k-7- Um]g . 


(A36) 


For vanishing mass m = 0, the mode functions 

1pk,s,e{x) = ■Ufc^s.e 6+“'" "’ 


(A37) 


are simply the limit m —)■ 0 of Eq. (A28), but the normal¬ 
ization of the spinor parts is changed: 


'^k.s — 


1 

71 


ifci 


6 


= 71 


\k\ (ELi 


(A38) 


This ensures that products of spinors are still finite. Ex¬ 
plicitly, 


^k,s ^p,r — 


1 - 


|fc||p| 


{2si{k^p^ — k^p^) + k ■ p) 


1 1 


2 |fc||p| 


[ 2 s(fc^p^ — k^p^) — i{k‘^p^ — k^p^) 5s,-r 


(A39) 


1 f2sp^ 2sk ^^ ^ 

2 

lfp^-2sip^ k^-2sik‘^\ s ^ An\ 

-2\^l - 


^k.s ^p^r — ^k,s ^p^r 


^k,s ^p,r (^fc,s ^p,r) 5 

and for k = p all products vanish: 

'^k.s.e ^k,r,6 — 0 • 


(A41) 

(A42) 

(A43) 


Note that these products are not simply the limits m —)■ 0 
of Eqs. (A32) and (A35). Finally, 


Yl ^isUk,s,B = ^ • 7b 


J2^tsVk,s,B = ^[^-7 


(A44) 


(1,1) dimension .— For finite mass m ^ 0, the mode 
functions can be decomposed as 


^k,s,s{x) = J Uk,s,e 6+“'"“ 


(A45) 


where the spinorial parts are 

lojk + m 


'^k.s — 




2m 


(s 


ujk +m 
2m 


cjfc+m ’®J 
fc 0-3^ 

uik+m'^ 

6 


(A46) 


Here, is the third Pauli matrix. The products of spinors 
are 


^p,r — 


^p,r — 


A/(wfc + m){ujp + m) 


2m 

X 

kp 

{ijJk + m){L 

C(wfc -1- m){LOp -1- m) 


(A47) 


2m 

2sp 2sk , _ / A 

X I - A- ^ - Ss,r A48 

ojp + m ojk + m 
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b. Elements of canonical quantization 

Canonical quantization of the free spinor field results 
in the mode expansion 

4'(t,£c) = '^ak^s'^k,s,+ {t,x) + bl^^'4)k,s-it,x) (A58) 

k.s 


in the Heisenberg picture, or 

4'(a:) = ^ak,s^k,s,+ {x) (A59) 

k,s 

in the Schrodinger picture. Here, ak,s and bk^s are par¬ 
ticle and antiparticle ladder operators. They satisfy the 
canonical anticommutation relations 

= Sk,pSs^r 1 5 (AGO) 

and all yet undetermined relations vanish. Again after 
renormalizing the vacuum energy, the Hamiltonian of the 
free spinor field is 

='^0Jk (al^^ak,s+ bl^^bk,s^ ■ (A61) 

k,s 

As before, particle number eigenstates are generated from 
the vacuum by application of creation operators. However, 
the order of the creation operators is now relevant, because 
they anticommute. 

Appendix B: On the normalization of field modes 

We briefly sketch why the mode functions of the scalar 
field are normalized with respect to Eq. (96), and those 
of the spinor field according to Eq. (97). 

The underlying reason is that the Lagrange density 
of the free scalar field <& given in Eq. (A3) contains a 
double time derivative of the field, which leads to the 
conjugate momentum H = The Lagrange density 

Eq. (A23) for spinor fields on the other hand, only 
contains a single time derivative of the field, and leads 
to the conjugate momentum H = ifkl. This in turn 
makes said normalization conditions necessary to ensure 
that the canonical (anti)commutation relations of the 
ladder operators are met. To make the connection evident, 
assume the modified relations 

[cik^ Up] — CX dk.p , ^ ^k,p ^s,r , (®1) 

and so on. 


1. Scalar field 

To see how the normalization arises naturally, quantize 
the scalar field directly in the Heisenberg picture. The 
task is to find a linear operator <I>(a;) which (a) solves 
the Klein-Gordon equation 

(□-I-m^)$(a;) = 0 (B2) 

as an operator equation, (b) satishes the usual commuta¬ 
tion relations with its canonical conjugate field H = i9o^i 
namely 

[$(a;),n(x)]|^^^^^ =iS{x-y) , 


(B3) 
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and (c) is self-adjoint: 

<^>\x) = <^{x) , tf(a;) = n(x). (B4) 

In the first step, assume that 4> can be expanded in a 
Fourier series using the mode functions Eq. (AlO): 

= X! ‘fikix) , n(x) = ^ nfc(t) tpkix) , 

k k 

(B5) 


The Hermitian conjugate of the field is automatically 
4'^(x) = — 1115 ( 0 :) by definition of 11. 

Now assume once again that the operator can be ex¬ 
panded in a Fourier series 


k 

115(0;) = ^ nfc,5(t) (pk{x), 

k 


(B15) 


where nfc(t) = do^k{t)- In the Fourier domain, the Klein- 
Gordon equation becomes 

(5o+Wfc)$fc(t) (B 6 ) 

and the canonical commutation relations imply 

[$fe(f),nt(t)] = i4,p , (B7) 

where we have already exploited that 

$l(t) = 4>_fc(t), ni{t) = n.k{t). (B 8 ) 


where 11^,5 (t) = 14^1^ 5 (i)- Translated to the Fourier 
domain, the problem becomes the following: find a linear 
operator 4'^(t) that solves the differential equation 

(z 7 % + 74,-m)^^(t) = 0 (B16) 

and satisfies 

{4/^(t),n_p,5(t)} = i4,p5^. (B17) 

In the next step, make an ansatz that again separates 
the functional part from the operator part: 


Solving the problem now amounts to finding 4>fc(t). 

The next step is to make an ansatz for 4>fe: 

$fc(t) = ukit) ak + u*_kit) a^_^. . (B9) 

The Uk{t) are simply functions, while Ofc are linear opera¬ 
tors. The Klein-Gordon equation reduces to 

(do + ^l) Ukit) = 0 , (BIO) 

which means that Uk is a linear combination of 
Requiring the commutation relations Eq. (Bl) for the 
ladder operators, the commutation relations Eq. (B7) for 
the field lead to the condition 

- ia [ukit)dou*k{t) - u*k{t)doUkit)] = 1 . (Bll) 


'^kit) = ^k,sit) Uk,s + v\^it) 

S 

nfc,5(i) = + U^-k,s,Bi^) 


(B18) 


The Dirac equation then reduces to 

{i-f°do + fh - m) Uk^^it) = 0 

(i 7 °ao - 7 *fc* - m) Vk^^it) = 0 . 


(B19) 


The spinors Eq. (A29) that solve the classical Dirac equa¬ 
tion satisfy 


(t^^m “ ua) Uk,a = 0 
+ m) Vk,s = 0 , 


(B20) 


In combination with {cpp , ipk) = Sk.p, we obtain 

{tpk j ‘Fp) — ^k,p j (®^^) 

a 

which, in case of the canonical commutation relations of 
the ladder operators a = 1 , is indeed the normalization 
condition Eq. (96). 


2. Spinor field 


so it is easy to see that the above equations are solved by 


Uk,si*) = CsUk,se 

UkA*) = -Ds e+‘“'=‘ , 


(B21) 


where Cs, Ds G C. Our ansatz thus reads 


'ftit) = 

Cs ak,s + e+-'=‘ . (B22) 


A very similar analysis for spinor fields leads to the 
normalization Eq. (97) for spinor mode functions. This 
time, the task is to find a linear operators 'I'(a;) which 
(a) solves the Dirac equation 

— m)'^{x) = Q (B13) 

as an operator equation, and (b) satisfies the canonical 
anticommutation relations with its conjugate field 11 = 

{4/^(x), n 5 (j/)}|^o^^„ = Ax -y)5i. (B14) 


The anticommutation relations Eq. (B17), together with 
those for the ladder operators Eq. (Bl), then lead to the 
condition 




Direct computation shows that this condition is met for 
Cs = Ds = -ym (awfc)”^. The field operator thus enjoys 
the mode expansion 


^'^(a:) = '^A,s,+ ix) ak,s + A,s-ix) b\ ,, 

k,s 


(B24) 
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with mode functions 

■ipk,s,+ {^) = 

These satisfy 


auJkL' 


auJkL' 




-ik-: 


, Xk^s e 


+i/i:-a: 


('0fc,s,e j ‘^p,r,s) — ^s,r ^e,(5 


(B25) 


(B26) 


and, choosing a = 1, we finally obtain Eq. (97). 


Appendix C: Proofs of Wick normal-ordering 
1. Quantized scalar field 

In this section, theorem VL2 is proven on the basis of 
theorem VI.l. We already know from theorem VI. 1 how 
to rewrite the time-ordered product of fields as a sum of 
normal-ordered products with more and more contrac¬ 
tions. To normal-order the entire product and obtain 
the expression stated in the theorem, we therefore merely 
need to commute all creation operators from right to the 
very left, and all annihilation operators from the left to 
the very right. In doing so, we will pick up commutators, 
between ladder operators and fields as well as between the 
different ladder operators themselves. By theorem VI.l, 
the full product 

{ap)p {bp)p [T<^>^^ (xi) • • • (x„) ] {al)k (4)fc 

can be reduced to a sum of products of the form 


the expression. Then, the creation operators have to be 
commuted past the field operators to the very left. Recall 
the commutation relations of the ladder operators of the 
quantized complex field Eq. (A 14), with all other pairings 
vanishing. In the first step, we need the resulting commu¬ 
tation relations between annihilation operators and the 
field parts. The four possible commutators are 

[ofc, 4>(x)] = 0 , [ttfe, $^(x)] = 1 ip*k{x) , 

[bk,<^{x)] = tiplix) , [ 6 fe,$'f(x)] = 0 . 

Consider the simplest case, when there is only one 
annihilation operator ak on the left, and no annihilation 
operators on the right. To move ak to the right, it has to 
be commuted with each and every field 'I''^*(xi) operator 
sooner or later—or more precisely with their constituent 
ladder operators. Each field operator decomposes into 
two part, one containing the annihilation operators, and 
the other the creation operators: 

ak :(xi) ■ ■ ■(xn) ■■ 

= ak : (xi)--- (x,) + (x*)] •••$"" (x„): . 

(C3) 

This decomposition can be done for all the other n — 1 
field operators as well. The result are 2”“^ summands 
containing as a factor (and all possible combinations 
of the parts of the other fields), and another set of 2 "“^ 
summands containing multiplied with the same set 
of combinations of parts of the other fields. Let us pick 
one such set and call it S for brevity. For example, S 
might be 

S' = (xi) • • • (xi_i)4>^‘+i (xi+i) • • • (x„) . 


(®p)p (^p)p ■ ^ ( 2 ^ 1 ) • • • $ " (x„): (al.)k (^^)fe • 

For these holds 

Lemma C.l. Using the notation of theorem VI. 2 , the 
following identity holds: 

{ap)p {bp)p ■ 4? ^ {xi) • • • $ " (x„): (aj,)fc {^k)k 

= :(op)p {bp)p (xi) • • • (x„) (a;[,)fc (&t )^ : 

I I 

+ 2^ ■ ap.^ ■ ■ ■ Op. • • • $ ^ (xj) ■ ■ ■ b^^ : 
single 

contractions 

+ (Cl) 

double 

contractions 

where eontraetions are only allowed between two ladder 
operators, and between a ladder operator and a field (not 
between two fields). 

Proof of Lemma C.l. To achieve normal-ordering, the an¬ 
nihilation operators have to be moved past the field oper¬ 
ators as well as the creation operators to the very right of 


Note that 

^S = :$^H 2 ;i)---$^)---$^"(cr„): , (C4) 

s 

where the term with hat is omitted. Since the entire 
product is normal-ordered at the end, and will 
end up at different positions inside S: 


:$''i(xi) •••$""(x„): 

= ^ [Si4/^-(x.)S2 + S*M^'{x.)S*, 
s 


(C5) 


where Si, S* are two differen partitions of S: S 1 S 2 = S 
and Si*S| = S. Thus, 


= • • • + afcSi$^- (x,)S 2 + akSl<I>l^' (x,)S* + • • • . (C 6 ) 

We now start moving ak to the right, and the contribu¬ 
tion from 4>(xi) will be the following: ak is commuted 
through Si (creating an additional term containing the 
commutator with other field parts in each step), and will 
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eventually end up next to Commuting and 

gives two terms 

Siak^'''{x^)S2 = [flfc, + Siak^'''{xi)S2 , (C7) 

where we have used that the commutator is proportional 
to the identity operator 1 in any case. The same is true 
for the second term containing : 

= K, (x,)]5 + {x,)S*2 . (C8) 

After combining the two commutators one obtains 

Ofc : ■ 

= -h [ofc, (a;i)]S' + ixi)akS2 

+ Sl<^l^'akS*2 + --- . (C9) 

All the terms containing commutators with the other 
fields have not been made explicit. Since this works for 
any combination S of the ladder operators of the other 
fields, the sum of all these will just give 

K, 4*"' (xj)]: (xi) • • • (xi) • • • $"" (x„): (CIO) 

from the first term according to Eq. (C4). The other 
two terms need to be processed further, giving rise to 
commutators with the other fields. In the end, after 
commuting oj, all the way through S there will be a 
commutator term like the above for every field, and a 
term which yields 


s 


ak 


= : 4>''i(xi) • • • $'"(x„) :afc 

= :afc$''i(xi) (x„): 


(Cll) 


according to Eq. (C5). By virtue of the commutators 
Eq. (C2) and the definition Eq. (162) of the contractions, 
we finally obtain 


Ofc : (xi) • • • 4>''" (x„): = : afe4>"i (xi) • • • (x„): 

I- 1 

+ ^ : afe4>"i(xi) • • • 4>"*(xi) • • • 4>‘'"(x„): . (C12) 

i=l 

Since it was nowhere used that a is a particle annihila¬ 
tion operator, the same relation holds for antiparticle 
annihilation operators b. 

Now, if there is a single creation operator on the right 
instead, the relevant commutation relations are: 

[^(a;), 4 ] = l‘/5fc(a;) , [$^(x),al,] = 0 , 

[$(x),4]=0, [$1’(x),5j[] = l(/?fc(x) , 

and a similar reasoning eventually gives 

: $"i(xi) • • • $""(x„) : = : ^''^(xi) • • • $''"(x„)a^ : 

” I- 1 

+ : $''i(xi) • • • $‘'’(xi) • • • $""(x„)a^ : (C14) 


in terms of the contractions Eq. (162). 

If both types are present, Eq. (C12) immediately allows 
to write: 


Ofc : • • • ^’""(x„): a], 

= : (xi) • • • (x„): 

" I- 1 

+ ^ : Uk^^^ (xi) • • • (xj) • • • $*'" (x„): a\, . (C15) 

i^l 

The first term on the right-hand side evaluates to 

:$"i(xi) (x„): alak'. 

+ :afc4>''Hxi) •••$''"(x„)4: (C16) 

using the contractions Eq. (161). The second term in 
Eq. (C15), as well as the first term in Eq. (C16) can then 
be evaluated using Eq. (C14). In total, one obtains 

Ofc : (xi) • • • (x„) ■■ al = : ak^^^ (xi) • • • 4*"" (x„)a^ : 

" I- 1 

+ ^ : afc4>"i (xi) • • • (xi) •••$"" (x„)a^ : 

+ ■■ ak^^^ixi) ■ ■ ■ 4>""(x„)4 : 

” I- 1 

+ X! • ( 3 ^ 1 ) •• • (a;i)4'"" (x„) • • • a],: 

n n I_I 

+T T • (a^i) • • • f"" i^j) • • • (xi) •••$"" (x„)4 : • 

i=i 3=1 ' ' 

(C17) 

Repeated application of Eqs. (C12) and (C14) finally 
yields Eq. (Cl) in the general case. 

□ 


Proof of Theorem VI.2. Combining theorem VI. 1 and 
lemma C.l directly allows one to express 

{<^p)p {bp)p [T4> ^(xi) • • • <I> "(x„)] (al,)k ( 4 )fe 

as normal-ordered, plus sums of the same normal-ordered 
products with all the contractions (single, double,.. .fully 
contracted), defined in Eqs. (155), (161) and (162). If the 
vacuum expectation value is taken, all terms which still 
contain uncontracted normal-ordered operators vanish. 
If the total number of operators (ladder operators and 
field operators) is odd, there will be at least one normal- 
ordered operator remaining and all terms vanish. Only 
if the number of operators is even, the terms where all 
operators are contracted survives since all the contractions 
are proportional to the identity operator. Note that some 
or all of these surviving terms can still be zero, since many 
of the contractions also vanish. □ 
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2. Quantized spinor field 

We now prove theorem VI.5. Analogous to the proof 
of theorem VI.2, we first prove the following 

Lemma C.2. Using the notation of theorem VI. 5, the 
following identity holds: 

(Op r^p^r i^p.r}p,r • 'I^ (^l) 

• • • 41 "■ (Xn) '■ i^k^g)k,s 

— •(npjr)p,r i^p.f^p.r 41 ^ (^l) ’ ' ’ 

■ ®(p.'’)i ''' ®(p.»’)i • • • 4* ^ (xj) ■ ■ ■ : 

single 

contractions 

+ E'" ■ (C18) 

double 

contractions 

where eontraetions are only allowed between two ladder 
operators, and between a ladder operator and a field (not 
between two fields). 

The proof is very similar to the proof of lemma C.l, 
but with the additional complexity that one has to keep 
track of signs introduced by anticommuting the involved 
operators. 

Proof of Lemma C.2. To achieve normal-ordering, the an¬ 
nihilation operators have to be moved past the field op¬ 
erators and the creation operators to the right, and the 
creation operators subsequently past the fields to the left. 
Recall the anticommutation relations between the ladder 
operators Eq. (AGO), with all remaining anticommutators 
vanishing. For the first step, we need to know how the 
annihilation operators anticommute with the fields. From 
the above relations and the separation of the spinor field 
according to Eqs. (173) and (174) one obtains: 

4'(a;)} = 0 , {afc,s, 4'(a;)} = 1 f’k.s.+ ix) , 
{&fc.s,4'(x)} = 0 , {&fc,s,4'(x)} = l'i/'fe,s.-(a^) • (C19) 

Let us begin with one annihilation operator on the 
left: to move it to the right, it has to be anticommuted 
with any field operator ’^'^'{xi) at some point, picking up 
an anticommutator in the process. Each field operator 
decomposes into two parts, one containing the creation 
operators, one the annihilation operators 

ak,s ■ 4'"^ (a:i) • • • 4'"' (x,) • • • 4'*'" {x„): 

= ak,s ■■ 4-*^^ (a;i) • • • [4'^* (x,) + (x,)] • • • 4''^" (x„): . 

(C20) 

The same can be done for the other fields, and we essen¬ 
tially end up with sums of products of ladder operators 
(with sums in front, and multiplied mode functions), where 
each term contains one ladder operator per field. Let us 


pick out one of these products of ladder operators from 
the field 4''^i (xi),..., 4'*^" (x„) without 4'*^’ (xi) and denote 
it S for brevity. Both parts of 4'‘^*(xi) are multiplied with 
S, but since the entire product of field operators is then 
normal-ordered, the position where they are inserted into 
S will be different. Through normal-ordering, 4'^‘(xi) will 
be shifted from position i to some position i + m, picking 
up a sign (—I)™, and 4'j’(xi) will be shifted from position 
i to some other position i + I, picking up a sign (—1)^ 
Therefore, 

Ofe.s : 4'‘'i(xi) • • • 4'"'(xi) • • • 4'""(x„): 

= ■■■+akA-^^Sl^'aHx^)S 2 

+ akA-^ysm'ix,)S^ + ... , (C21) 

where 5'iS'2 = S = 8 ) 82 . Consider the first term, con¬ 
taining 4'^*(xi): the ladder operator ak,s is successively 
anticommuted through the operators of ^i, creating an ad¬ 
ditional term containing the anticommutator, and picking 
up a sign for ak,s in each step. By the time ak,s reaches 
4'^‘(xi), it already had to be anticommuted through the 
(i + m — 1) normal-ordered ladder operators in 81 and 
has therefore picked up a sign (—1)*+™“^, which leaves 
us with an overall sign (—l)*+2 "‘-i = Anticom¬ 

muting ak,s and 4'^*(xi), we get two terms: 

i-iy-Aiak,snA^)S2 

= Aiy-Aak,s,nAy}8- 

{-iy-^ 8 i'Sl'{xi)ak,s 82 , (C22) 

where we have used the fact that the anticommutator is 
proportional to the identity operator in any case, and can 
therefore be pulled in front of the product S'! 5*2 = 8. The 
same argumentation works on the second term, 4'^’(xi), 
and we obtain 

{-iy-AAk,sn'{xy 8*2 

= i-iy-Aak,s,nyxy}8 

- {-iy-Amyxi)ak,s8A (C23) 

The two anticommutators can be combined to 

{ak,sAyixi)} + {ak,s,%yxi)} = {afe,s,4'"*(xi)} (C24) 

such that Eq. (C21) becomes 

ak,s ■■ 4'‘'^(xi) • • • 4'"‘(xi) • • • 4'""(x„): 

= • • • + (-l)*5i4'^*(x,)afc,,A2 + {-iy 8 l%'{xyak,s 8 ^ 
+ (-l)*"'{afe.„4/^‘(xi)}5+... . (C25) 

We have not written out the anticommutators ak,s has 
picked up with the operators in and 8 ). Since this 
needs to be done for any of the field operators, we end 
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up with 

ak,s ■ : 

= (-1)” : • • • 4'*'" (a;„): Ofe,* 

n 

+ : ^'"1 {xi) ■ ■ ■ 

i=l 

• (C26) 

The first term on the right-hand side is normal-ordered, 
and the sign can be absorbed into anticommuting the 
annihilation operator back to the beginning of the product. 
The second term can be rewritten by taking into account 
the anticommutation relations Eq. (C19) and defining 
contractions like in Eq. (186) 

ak,s ■ (a;i) • • • 4'"" (a;„): = : (a^i) • • • 4'*'" (xn) ■ 

I I 

i=l 

Note in particular that the sign is absorbed into the 
definition of the contraction. Since we have nowhere 
used that a is a particle annihilation operator, the same 
identity holds for antiparticle annihilation operators b. 

Now consider what happens if a creation operator is 
added at the right. Using the last two equations, we can 
immediately write 

ak,s ■ 4'"^ (xi) • • • 4'"" (xn) : a], 

= (- 1 )” : 4'"! (xi) • • • 4'"" (xn ): afe.sflp^r- 

” I- 1 

+ X! :afc,s4'"i(a;i) • • • 4'''*(a;i)a],_,,4'''"(a;„): . 

(C28) 


Taking into account the anticommutation relations be¬ 
tween the ladder operators and defining contractions 
Eq. (185), the first term on the right-hand side can be 
expressed as 

: ak,s'^''Hxi) ■ ■ ■ 4'""(x„)d^ ,.: 

+ (-l)"+i: T^i(xi) + • • • T^’*(cr„): al,ak,s ■ (C29) 

To handle the remaining terms, we need to evaluate ex¬ 
pressions of the form 

: 4/^1 (xi) • • • 4/^* (x,) • • • 4/^" (cr„): . (C30) 

The relevant anticommutation relations are 

{'^ix),al j = 0 , {4'(a;),4 J = lV'fc.s.-H(a;) , 

. _ . _ 

= 0 , {4'(a;),6l_J = l'ipk^,_{x) , 

and a reasoning analogous to the one that led to Eq. (C27) 
eventually yields 

: 4'"i(a;i) • • • 4'""(a;„): a^ ,. 

= : 4'"i(xi) • • • 4'""(a;„)a^ ,.: 

” I- 1 

+ J2-'^"Hxi)--- (xi)-^^- (x„) • • • al, : . (C32) 

i^l 

Using this in Eqs. (C28) and (C29), one obtains 


ak,s : 4'"i(a;i) • • • 4'""(a;„): = : afc,s4'"i (a;i) • • • 4'*'" (a;„)a^:-f-^ : afc,s4'‘'i (a;i) • • • 4''"* (cci) • • • 4'"" (a;„)a],: 


2 = 1 


+ : afe,s4'"" (a^i) • • • 4'*'" (a;„)a],_,,: + E ’ (a^i)''' 4'"' (a;*)4'"" (x„) • • • a],: 


2=1 


+EE : afc,s4'"i (xi) • • • 4'"^ {xj) • • • (xi) • • • T"" (x„)a],,, 


(C33) 


i=ij=i 


Iterative application of Eqs. (C27) and (C32) then 
yields Eq. (C18) in the general case. □ 

Proof of Theorem VI. 5. Similar to the proof of theo¬ 
rem VI.2, combining theorem VI.4 with lemma C.2 yields 
theorem VI.5: at first, one obtains a completely normal 


ordered term plus all normal-ordered terms with single 
contractions (between two fields, between two ladder op¬ 
erators and between ladder operator and field) plus all 
normal-ordered term with two contractions and so on. 
Taking the vacuum expectation value, all terms vanish 
due to the normal-ordering, with the exception of those 
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terms that are fully contracted and therefore proportional 
to the identity 1. Full contraction is only possible for an 
even number of operators, such that products containing 
an odd number of operators vanish. □ 


3. UDW-type detector 


Theorem VI.7 is essentially a simplification of the the¬ 
orem VI.5 proven in the previous section. Again, first 
consider 


Lemma C.3. Using the notation of theorem VI. 7, the 


following identity holds: 

= ■■[cr~r■ ■ ■ fi{tn)[a+]’’■■ 

+ ^ ; (T~ • • ■ fi{ti) ■ ■■ il{U) ■ ■ ■ fj,(t„) • • • cr+ : 

single 

contractions 

+ ^ : ct" • ■■ fifti) ■ ■ ■fJ.itj) ■ ■ ■ li{ti) ■ ■ ■ g.{t„) • • • cr+ : . 

double 

contractions 

(C34) 

Proof of Lemma C. 3. The proof runs analogous to the 
proof of lemma C.2, but is somewhat simpler since there 
is only one type of ladder operator (ct) instead of two 
(a, b), and the monopole operator is self-adjoint, unlike 
the spinor field. This time, the relevant anticommutation 
relations are Eq. (124), and 

{a-,/r(t)} = a+} = 1 (C35) 

instead of Eq. (C19) and Eq. (C31). This gives rise to the 
contractions defined in Eqs. (198) and (199). Beyond that, 
the proof is follows the exact same steps as before. □ 

Proof of Theorem VI. 7. The theorem follows directly 
from theorem VI.6 and lemma C.3, analogous to the 
proof of theorem VI.5. □ 
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